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connected financial systems. However, beyond contagion research, these topics are 
also essential for risk-neutral pricing. In this article, we therefore study and com¬ 
pare valuation algorithms in the standard model of debt and equity cross-ownership 
which has crystallized in the work of several authors over the past one and a half 
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1. Introduction 


Since the turn of the millennium, research interest in systemic financial risk has steadily grown 
with a noticeable pick-up in the number of publications ov er the past five y e ars. One main field 
of inte r est is defau l t cont agion , see for ins t ance t he works of Acemoglu et ah ( 20131 1. Elliott et al 


1 2ni3l l . Gai et al. (2011) and Nier et al. (2007), or Staurn (12 012 ) for a survey. Ma ny publica- 


were 


tions regarding systemic risk refer to the seminal work of lEisenberg and ( 2001 ). who 
the first to structurally model financial systems in which firms can hold each other’s financial 
obligations as assets under the as sumption of limited liability. However, the core idea behind 
the mo d el of lEisenberg and Nod , which can be interpreted as a multi-firm extension of the 
MertonI (jl974l ) model where cross-holdings of zero-coupon debt between members of the system 
is allowed, received somewhat less attention by the research community. The main difference 
between Eisenberg and Noe and the standard multi-firm Merton model is that prices at maturity 
are not trivially determined since the value of one firm’s equity or debt m ay depend on the value 
of the debt of any other firm in the system. lEisenberg and Nod (j200ll ) gave conditions under 
which only one equilibrium solution exists at maturity. Together with a finite numerical algo¬ 
rithm that they provided, the model could not only be used for default contagion rese arch, but 
also fo r risk-neutral n o-arbitrage valuation under financial interconnectedness (see also iFisched . 
20141 1. lElsingeu (j2009j) generalized the Eisenberg and Noe setup by also including cross-holdings 
in equity, and by allowing a seniority structure of the liabilities. A numerical algorithm was pro¬ 
vided, however, a finite number of steps to the equilibrium price vector could not be guaranteed 
anymore. 
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Already in 20n3. Teruyoshi Suzuki had - unbeknown to him - generalized the 


Eisenberg and Noel setup to the situation where debt of one si ngle seniority and equity could 


be cross-owned within a financial system (Suzuki, 2002). Unlike Eisenberg and ( 2001 ) and 
Elsinger ( 2009l b Suzuki , who had the clear intention of generalizing Merton ( 1974l l. provided 
a Picard Iteration as the numerical means of calcu l ating price equilibria. In a further gen¬ 
eralization of Suzuki’s and Elsinger’s work, Fischer ( 2014l i extended the structural model of 
interconnectedness to the case where liabilities cou ld be de r ivativ es in the sense of a dependence 
on other liabilities or equities in the system. As in ISuzukil ()2002l h the numerical procedure pro¬ 
vided to solve the liquidation value equations at maturity was the Picard Iteration. However, 
unfortunately, a Picard Iteration cannot warrant an exact solution in finitely many steps. 

So, while there exists a small but growing amount of research on the existence and the unique¬ 
ness of price equilibria in systems with financial interconnectedness, the provided algorithms 
mainly reflect the individu al authors’ particular ap proach to the problem. For instance, there 
also exists a publication bv iGourieroux et al.l ( 2013l l which mentions that, in the Elsinger model 
with two debt seniorities but no equity cross-holdings, a simplex method can be applied. How¬ 
ever, comparative studies of the different methods seems to be absent from the existing literature. 
Furthermore, at present, n o numer i cal algorithm for the setup with cross-holding s of equity and 
one seniority class of debt ( Suzuki . 2002 : Elsinger, 20091 : Gourieroux et ah . 2012) is known that 
reaches the exact solution in a finite number of calculation steps. 

For these reasons, the article at hand has three main objectives. First, we want to provide an 
overview of the already existing valuation algorithms by unifying notation and by embedding 
them in one general model framework. Second, we provide a new type of algorithm which is a 
hybrid of Eisenberg and Noe’s and Elsinger’s approach that has improved convergence properties. 
Third, we introduce a whole range of algorithm versions which are based on the three different 
types - namely Picard, Elsinger, and Hybrid - which will reach the exact solution (if existing 
and unique) in a finite number of calculation steps. Introducing these new algorithms, we 
show that for the three types of algorithms there always exists an increasing and a decreasing 
version - depending on a properly chosen (and explicitly given) starting point. Furthermore, 
we use default set techniques and linearization techniques to achieve finiteness. A simulation 
study hnally allows to draw some conclusions about the efficiency of the presented numerical 
algorithms with respect to model parameters such as system size. 

The structure of this paper is as follows. In the next section, we will establish the model 
and necess ary assumptio ns fo r a uniqu e solut ion of the financial system. The existing valuation 
methods of Suzuki ( 2002) and Elsinger (1200911 a re pr esented in Sectio n [3l where a hybrid version 
of the algorithms of Eisenberg and Noa ~( 200ll i and Elsinger ( 2009l i is developed as well. The 
algorithms of this section are all non-finite. In the fourth section, we introduce a new class 
of valuation algorithms based on default set techniques that hnd solutions in finite time. A 
simulation study in Section 0 compares the runtimes of the different algorithms for different 
classes of financial systems. In Section [6l we conclude. A technical appendix follows. 


2. Notation and Model Assumptions 

For two matrices M = (Mjj)jj=i^,,,^„ € and N = E R”-^" we write M > N 

if Mij > Nij for all i,j E {1,... ,n} and M > N if Mij > Nij for at least one pair For 

two vectors u = (ui,... ,UnY E R'^ and v = (ui,... ,VnY E R” the definition of u > v and 
u > V is analogous to the conventions for matrices above. A matrix M E is said to be left 

substochastic if Mij > 0 for all z, j E {1,..., re} and if YYd=i ^ij — ^ J S {1,..., re}. The 

symbol I„ is used for the re x re-identity matrix and 0^ is used for a (column) vector of length re 


3 








































































that contains only zeros. Additionally, O^xn stands for an (n x n)-matrix with only zero entries. 
For a vector u € M"", the expression diag(u < 0,i) stands for an (n x n)-diagonal matrix where 
the i-th entry is 1 if Uj < 0 and 0 otherwise, i.e. 


diag(u < 0 „) 


1, for i = j and Ui < 0, 
0, else. 


( 1 ) 


All operations such as the minimum, min{-}, the maximum, max{-}, or the positive part (■)'*' 
are applied element-wise to vectors and matrices. The norm in this paper is the £^-norm on M"" 
defined as 

n 

||x|| := ||x||i = \xi\ for x € M"". (2) 

i=l 

The corresponding norm for a left substochastic matrix M £ is given by 

n 

||M|| := ||M|| 1 = max ||Mx|| 1 — max ^ ^ ]\^ij ^ 1, (^} 

ii^ii 1 ^ 


meaning that ||M|| is the maximum of the column sums. One easily can show that ||Mx|| < 
||M||||x||. 

We consider a system of n financial entities, and denote M = {1,... ,n}. In the following 
these entities are simply called “firms”. Each firm owns exogenous assets, that are defined in 
the next step. 


DEFINITION 1. Let Oj > 0 denote the market value of the exogenous assets held by firm i. 
As the name implies, these assets are priced outside the considered system in the sense that the 
eapital structure of the n firms has no influence on the prieing mechanism of such an asset. By 
a = (oi,..., OnY € we denote the (eolumn) veetor of the exogenous assets. 


Moreover, we assume that the firms have outstanding liabilities with a nominal value at 
maturity of di for each firm i. These liabilities are summarized in the vector d £ (Mg')"'. In our 
framework we assume that the entries of d are constant. Since it is assumed that the exogenous 
assets’ prices are given by the constant vector a, the results of this paper presented in the 
Sections [2] toll also hold if d depends on a, i.e. if d = d(a). However, for the remainder, we 
will write d for convenience. This definition of the liability vector allows the interpretation that 
the di are simple loans or zero coupon bonds since they are not derivatives that can depend 
on the other assets within the system. T he case of con s tant liabilities is u s ed in most existing 
public ations in this field, see for example Suzuki ( 2002l l. Gourieroux et al. ( 2012 1 and Elsingei 
( 20091 1. The more general ca se in which d d epends for example on the endogenous assets, is also 
treated in the literature (see Fischer . 2014l l. 

To take the interconnectedness of the firms into account, we allow that each firm can own a 
fraction of the liabilities of the other firms. To formalize these possible cross-holdings, we use 
ownership matrices. 


DEFINITION 2. The left substochastic matrix in which the entry 0 < M'j < 1 

denotes the fraction that firm i owns of the liability of firm j is called debt ownership matrix. 
Since no firm is allowed to hold liabilities against themselves, we assume = 0 for all i £ Af. 
The entries Mfj of the (left substochastic) equity ownership matrix M® £ are defined as 

the fraction that firm i owns of firm j’s equity. 
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Note that the diagonal entries of M® must not be zero; that means it is allowed that firm 
i holds its own shares in which ca s e Mf- > 0. For the debt ownership matrix it is a common 
convention (cf. Eisenberg and Noe ( 2001 ) or Elsingei ( 2009l ll that = 0 since a firm cannot 
have debt obligations to itself. The tuple J- = (a, d, M®) is in the following sometimes 
referred to as the financial system. 

Associated with the liability vector, we consider the recovery claim vector r G The 

recovery claim vector represents the actual payments of the firms at maturity, i.e. in general 
we have since default risk is present. The value of the debt claim that firm i has 

against firm J is hence given by Mf-Vj and the total value of firm i’s debt claim against the other 

Lj J 

members of the system is Furthermore, denote by s G (IRq')"' the equity values of 

the n firms which means that the total recovery value of all system-endogenous assets that firm 
i owns is given by the i-th entry of 

M'^r -h M®s. (4) 

The basic assumption for the model is that equity is considered to be the residual claim 
which means that any outstanding liability has to be paid off completely before the shareholders 
receive a positive payment. Hence, the equity value s* of firm i is positive if and only if firm i can 
fully satisfy all their obligees. The Absolute Priority Rule immediatel y leads to the following 
liquidation value equations for the recovery claims and the equities (cf. Fischer . 2014l l: 


' = min{d, a -|- -|- M®s} 

. = (a-hM‘^r-hM®s-d)+, 


(5) 

( 6 ) 


where the sum in ([5]) contains no (•)+ since Theorem [T] will show that all solutions of this system 
are non-negative. A solution for the liquidation value equations ([5]) and ([6]) is hence a fixed 
point of the mapping <1> : where R = (r*,s*)* G (Rq")^” and 


4.(R) = 4. Q 


/ min{d, a -b M'^r -b M®s}\ 
V(a-bM‘ir-bM®s-d)+ ) ' 


(7) 


We will sometimes refer to the debt component of R and mean in such cases the first n compo¬ 
nents of R that represent the debt payments of the systems. The components n -b 1 to 2n of R 
we also call equity component for the same reasons. We are interested in finding the fixed points 
of <h, which we will also call solutions of the financial system iF = (a, d, M®). Without 
further constraints it is possible that there exist several fixed points. To ensure that the solution 
is unique, we have to make an additional assumption in which we need another property of an 
ownership matrix. 

DEFINITION 3. An ownership matrix M G R"'^” possesses the Elsinger Property if there 
exists no subset J F N such that 


Mij = 1 for all j G J. (8) 

i&J 


The name of this property is chosen because lElsingeii ([20091) is, by the best knowledge of 
the authors, the hrst one to use this assumption in the context of ownership matrices and the 
valuation of systemic risk. For our model, we demand that the considered ownership matrices 
fulfill this property. 


ASSUMPTION 1. The Elsinger Property holds for the debt and the equity ownership matrices 
and M® . 
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Note that the fact that and M® are hol ding mat r ices is equivalent with the existence 
of (I„ — and (I„ — as shown by Elsinger ( 2009l i. Moreover, Theorem [1] below 

will show that Assumption [T] ensures that there is only one fixed point of <h. To show this, we 
introduce the two vectors 


R 


great 


_ / fgreat \ _ 


d 


Weaty V(In-M®)-l(a + Mdd-d)- 


(9) 


and 


R, 


'Small 


f I'small 
V^small 


min{d, a}\ 

. (a-d)+y 


( 10 ) 


The vector Rgreat assumes that the debt payments are fully recovered so that in the debt 
component r = d. Note that even if for a fixed point R* = (g*) of it holds that r* = d, 
it must not necessarily hold that Rgreat = R*- The second vector Rsmaii emerges when the 
liquidation equations m and ([6]) are applied and the ownership structure of liabilities and 
equities is completely ignored. In this case the term M'^r + M®s that represents the income of 
each hrm stemming from debt and equity cross-ownership is set to zero. Hence, the firms only 
have the exogenous assets a as an income. The vector Rsmaii results from applying the mapping 
<I> in equation ([7D to the vector 02 n, i-e. 


LEMMA 1. With the definitions above, it holds that ( [Rsmaiij Rgreat]) C [Rsmaib Rgreat]• 

Proof. Assume R € [RsmaibRgreat]- Because of ifTTI) and the monotony of $, we have that 
$(R) > Rsmaib By definition of and Rgreat) only the lower n lines of the vector inequality 
ch(R) < Rgreat need to be shown. By Lemma IA3I and ([9l) , it holds that 

Sgreat “ M®Sgreat = (a M'^d - d)+ (12) 

and therefore 

Sgreat > (a + M'^d + M®Sgreat “ <i)^- (13) 

□ 


Before showing the importance of Rgreat and Rsmaii as upper and lower bounds of the solution 
R*, we need to introduce the terms default set and default matrix. For r > 0^ and s > 0^ the 
set 

(14) 

is called default set under r and s because - given r and s - the firms in D{r, s) are not able to 
fully satisfy their obligations and hence are in default. We say that hrm i is in default under 
r and s if i G Zl(r, s). For R = (r*,s*)* we will sometimes abbreviate the default set as D(R). 
The default matrix corresponding to r and s, A(r,s) G is dehned as 

A(r, s) = diag(a -|- M'^r -|- M®s — d < 0„) (15) 

and is the diagonal matrix with entry 1 for hrms in default under r and s at the corresponding 
position and with the value 0 for hrms not in default. With the new notation, we can show the 
crucial limiting property of Rgreat and Rsmaii- 


r n 

D{r, s) = < i G A/”: Oi ^ 
[ i=i 
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PROPOSITION 1. Let R* be a non-negative solution of the fixed point problem defined in 
dH). Then R* G [Rsmall, R-great]- 

Proof. Because of (fTT|) and the monotony of R* > Rsmall) so we only show the validity of 
the upper bound Rgreat- Since R* is a fixed point of we can write 





/ min{d, a + M'^r* + M®s*}\ 

V (a + - d)+ ) 



(16) 


Obviously, r* < d = rgreat, hence we reduce our considerations to the equity component of R* 
which, together with A(r*,s*) = A*, can be presented as 

s* = (a + + M^s* - d)+ = (I„ - A*)(a + + M^s* - d). (17) 


Because of (I„ — A*)s* = s* we can reformulate CZl) into 

s* = (R - A*)M®s* + {In - A*)(a + - d) 

= (I„ - A*)M^(I„ - A*)s* + {In - A*)(a + - d). 

Rearranging yields to 

s* = {In - {In - A*)M^(I„ - A*))-Hl„ - A*)(a + M\* - d). 

Together with Lemma lA4l in the Appendix, this leads to 

s* = (I„ - {In - A*)M^(I„ - A*))-nin - A*)(a + M^^r* - d) 

< (I„ - {In - A*)M^(I„ - A*))-1(I„ - A*)(a + M^d - d) 

< {In - (In - A*)M^(In - A*))-\ln - A*)(a + M^^d - d) + 

< (In - A*)(In - M^)-i(In - A*)(a + M^^d - d)+ 

< (In - M^)-\a + M'^d - d)+ 

— Sgreat? 

from which the assertion follows. 


(18) 


(19) 


( 20 ) 


□ 


Using the results from above, we can now show that there is only one fixed point which is in 
the interval [Rsmall, Rgreat]- The proof of the following theorem is given in the Appendix. 

THEOREM 1. Under Assumption{J\ and for an arbitrary financial system J- = (a, d, M®), 
there exists a unique fixed point of the mapping <i>. The fixed point R* is non-negative and 
R € [Rsmall, Rgreat ] • 


In the sequel, we assume that Assumption[T]holds so that J- has only one solution R* G 


i.e. 


R* = 


= $ ( = $(R*). 


\2n 


( 21 ) 


The requirements of Assumption [T] ar e less st r ict th a n the as s umpt ion that both ||M^|| <1 an d 
||M®|| < 1 that is used for example in iFischeil (j2014l l , Suzuki (j2002l l or lOourieroux et aP (|2012l l . 
This follows by the fact that in case of having the Elsinger Property, it must not necessarily 
hold that IlM'^ll < 1 and hence is no strict contraction anymore. However, the assumption 
still guarantees that the solution of the system is unique. 


7 
















3. Non-Finite Algorithms 


In this section, two existing solution algorithms that can be found in the literature are presented. 
One algorithm consists of the iterative use of the mapping <1> on a chosen starting vector and 
is given i n the first subsection. A modification of this Picard Iteration is used in the work of 
Elsingeu ([20091), where for the determination of the equity component, a more sophisticated 


subalgorithm is used (Section 13 21). I n the last subsection, a new al gorithm is developed that 
combines the ideas of Elsingeif l 2009l i and Eisenberg and Noe (2001) which results in a faster 
convergence of the procedure. 


3.1. The Picard Algorithm 

The most intuitive way to calculate R* for the system T consists of the iterative use of <h. It 
will be shown in this section that with an arbitrary starting vector G (M(j')^”', 

R* = lim $'(R°) = lim $ o ... o $ (R°), (22) 

i 


which is commonly known as the Picard Iteration. Since R* > 0, the range for the starting vector 
R° can be reduced to only non-negative vectors. Beyond that, the search for an optimal starting 
point can be limited to the interval [Rsmaih Rgreat], as shown in Theorem[T] A direct consequence 
is that any iteration procedure that aims to calculate R* should make sure that (i) no starting 
point of the iteration is chosen outside the interval [Rsmaih Rgreat] and that (ii) every interim 
result of the procedure also needs to be in that interval. Otherwise, the procedure is inefficient. 
For these reasons, we present an algorithm that can start with both, Rgreat and Rgmaii- 

ALGORITHM 1 (Picard Algorithm). 

1. For k = 0, choose R*^ G [Rsmaih R-great] CLnd e > 0. 

2. For k>l, determine R^ = ^>(R^“^). 

3. //llR'^-i - R*^|| < s, stop the algorithm. Else, set k = k + 1 and proceed with step [M 


We will use t he two expre ssion s Picard Itera tion and Picard Algorithm synonymously for 
Algorithm [TJ In Suzuki ( 2002 1 and Fischer ( 2014l l the Picard Iteration is the algorithm of choice 
to determine solutions of ([5]) and (|6|). 


PROPOSITION 2. In case of HP = Rsmaih Algorithmic generates a sequence of increasing 
vectors R^, and for R*’ = Rgreat o sequence of decreasing vectors. For all starting points, the 
algorithm converges to the solution R*. 


Proof. Let R*^ = Rsmaih then 


\ - f mm{d, a -h M'^rsmall + M®Ssmall}\ ^ f min{d, a}\ _ 
(Rsmallj ^ M^rsmall + M^Ssmall - d)+ y) “ ^ (a - d)^ 


'Small • 


(23) 


From the monotonicity of <h, it follows that for all iterates we have > Il^,k > 1. For 























= Rgrcat) first check that because of Sgj-eat = (In — M®) ^(a + M'^d — d)"'' and Fgreat = d, 


(a + IVI Tgreat “k hd Sgi-gat d) 

= fa + M'^d - d + M®s 


great 


— S, 


great 


+ S, 


great 


— f a. + ]VI d. d (1^2 ]V[ )Sgreat H“ Sgreat 


(24) 


= a + M'^d - d - (a + M'^d - d)+ + s 


great 


<o„ 


^ (Sgreat)^ — 


>great 


and thus 


<I>(Rj 


_ f min{d, a + M'^rgreat + M®rgreat}^ . 
V (a + M'^rgreat + M®Sg,eat - d)+ “ 


d 

Sgreat 


= R 


great • 


(25) 


Again it holds, due to the monotonicity of <h, that R*^+^ < R^,A: > 1. Hence for any R G 
[Rsmaih Rgreat] h follows because of R < Rgreat that $(R) ^ ^(Rgreat) < Rgreat cind with the 
same argumentation it follows that 4’(R) > ‘h(Rsmaii) > Usmaii- This means that any series 
from the Picard Iteration with a starting point in the interval [Rsmaib H-great] is bounded from 
above and from below. Since 4> is continuous, it follows that the series must converge to some 
R such that <I>(R) = R. According to Theorem [U there is only one fixed point, so it must hold 
that R = R*. □ 

Regarding the P icard It e ration wit h the starting point s Rsrnaii or Rp-rsar. it should be men¬ 
tioned that besides Suzuki ( 2002 1 and Fischer ( 2014l l. also Shin ( 2006l l considers a Picard itera¬ 
tion in a system valuation context. Shin’s model is one with debt cross-ownership and multiple 
seniorities, while eq uity cross-ownership is not consider ed. As su c h, the model is situated some¬ 
where between the Eisenberg and ( 200ll l and the Elsingei ( 2009l l framework. Instead of 
maturity values (as done here). Shin directly considers risk-neutral values at time 0, and takes 
for the start of the iteration procedure either a “conservative” viewpoint, where debt is assumed 
to have the value zero, or an “optimistic” viewpoint, where the value of debt is assumed to be 
the face value. Shin’s starting points therefore seem to be risk-neutral time zero equivalents to 
the here considered Rsmaii and Rgreat- 

The Picard Iteration - or any other iterative algorithm in this section - might not reach the 
solution R* in finitely many iteration steps. Examples of financial systems with this property 
can easily be constructed. Prom a computational or practical point of view this means that iter¬ 
ative algorithms like the Picard Iteration have the disadvantage that under some circumstances 
many iterations are needed to approach to R* sufficiently close, which makes these algorithms 
somewhat inefficient. The Trial-and-Error Algorithms presented in Section 0] do not have this 
drawback since for these procedures it is assured that they will reach the solution in a finite 
number of steps. 


3.2. The [Elsingeil Algorithm 

InElsmger ( 200911 . an algorithm for R* is presented which differs from the Picard Iteration. This 
procedure consists of splitting the two components of R, the equity and the debt component, 
and apply different computation methods on both components in each iteration step. For the 
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equity component, a sub-algorithm is applied where the equity payments of the system are 
determined assuming a fixed amount of debt payments. Denote this vector of debt payments in 
the following by f, hence 0,^ < f < d. Aim of the sub-algorithm is to find a fixed point of the 


mapping : 


with 


^>®(s; r) = {a + M'^f M®s - d) 


(26) 


This mapping represents the equity component of <1> for a given debt payment of r. The fixed 
point of <l)®(-;f) is denoted by s(f), i.e. 


<h®(s(r); r) = (a -|- M'^r -|- M®s(r) — d)’*' = s(r). 


(27) 


As shown in Elsingei ( 20091 1 . this fixed point exists and is unique since M® has the Elsinger 
Property. 

The following algorithm delivers for given f a series of vectors € M"" that converge to a 
vector whose positive part is the fixed point of (1261) . To explain this in more detail, first define 
for a given vector w € M"' the set 


P(w) = {i G : Wi > 0} 


(28) 


and the matrix 


r(w) = diag(w > Oj; 


(29) 


as the corresponding diagonal matrix. Note t hat these definitions of T’(w) and r(w) slightly 
differ from the original ones in Elsinger ( 2009l l. where a strictly larger sign was used. By our 
definition of default in (|14p . a firm with zero equity value can still be not in default in the sense 
that all obligations can fully served. This situation is referred to as borderline firms (cf. Section 
However, this modification does not change the forthcoming theoretical results. 


ALGORITHM 2A. 

1. For k = 0, set = a -|- — d and determine P(w^) and r(w^). 

2. For k>l, solve 'l'^fc(w) = w where 

T.^^fc(w) = w° -|- M®r(w^)w (30) 

and denote the solution by i.e. (w^"*"^) = Determine P(w^'*'^) and 

r(w^+i). 

3. If P{w^) = stop the algorithm. Else, set k = k + 1 and proceed with step\^ 

Before the properties of Algorithm I2AI are shown, we give some explanations for a better 
understanding of its functioning. The starting point is which is the difference between 
a + M^f and d. The sum represents the firms incomes on their balance sheet that consists of 
the external assets and the payments due to cross-ownership of debt. Note that in this step the 
potential income from equity cross-ownership is ignored since M® does not appear. The idea is 
now as follows; The firms not in P(w*^) are not able to fully satisfy their liabilities (assuming 
debt payments of f) and will be in default. On the other hand, the firms that are in P(w^) will 
be able to satisfy their obligees and can be regarded as solvent (again assuming debt payments of 
f), even though no intersystem payments due to equity cross-ownership are taken into account. 
As a consequence, the equity payments of the non-defaulting firms are added into the system 
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via the product M®r(w^)w. We can interpret the vector w^, as well as the other iterates w^, 
as pseudo equity vectors that give us information about solvent and defaulting firms under the 
current debt and equity payments. The fact that the entries of can be negative prevents 
that they can be naturally interpreted as equity vectors which is why we use the term “pseudo”. 

The difference compared to the Picard Algorithm is that a linear equation system is solved to 
achieve a new equity payment vector instead of applying to (f*, This is because for 

the fixed point of it holds together with (f30]l that 

= (I„ - M®r(w*^))"V° (31) 

Note that the inverse matrix exists since M® and hence M®r(s^) have the Elsinger Property. 

The vector can be interpreted as an “updated” version of since the equity of the 
non-defaulting firms that are in P(w*^) is included in w^. Based on the updated vector it 
might appear that some firms that are not in P(w*^) have now non-negative entries in w^. This 
can be concluded from that we will show later. But these firms are now also able to 

contribute equity payments to the system. Consequently, the system has to be updated again 
by determining w^. The procedure continues until the set of defaulting firms stays the same 
from one iteration step to the next one. 

PROPOSITION 3. Given a fixed vector of debt payments f > 0,i; 

(i) Alaorithm \2A\ generates an inereasing sequenee of vectors w^. 


(ii) Let 1 < I < n such that 

I := min{j € {0,1,..., n} : P(w^) = 

Then s(f) = is the fixed point of the mapping d>®(-;f). 


(32) 


(in) Let do = |P(w*^)| E {0,1,... ,n} be the number of firms with a positive entry in The 
fixed point s(f) is reaehed after no more than n — do iteration steps. 

Proof. (i) This part of the Proposition is shown by Elsinger ( 20091 1. We give a different 
version of the proof. Because of (1311) . the fact that r(w°)w° > 0,1 and using the series 
representation of (!„ — M®r(w®))“^ as shown in Lemma lA3l of the Appendix we get 


(33) 


= (I„ — M®r(w'^)) 

= (i„ -f M®r(w°) -F (M®r(w°))2 -h.. .)w° 

= w° -KM® r(w°)w° -MV[®r(w°)M® r(^v°)w° -f ... 


> 0 „ 


> 0 „ 


> W°, 


which is the induction start. Eor the induction step we assume > ^ and r(w^) > 

r(w^“^) following from it. We need to show that > w^, or, equivalently, = 

-I- e where e > 0^. Since r(w^)w^ > r(w^“^)w^ and -|- M®r(w^“^)w*^, it 

follows that 

u := -I- M®r(w^)w^ — = M®(r(w^) — r(w^“^))w^ > 0„. (34) 

With this definition we have that 

-b e = w° -h M®r(w^)(w'' -b e) = w° -b M®r(w*=)w'' -b M®r(w'=)e (35) 
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and we can rearrange to 

e — M®r(w^)e = + M®r(w^)w^ — = u > 0^. (36) 

Solving this for e leads to 


e = (I, - M^r(w^))-^u > 0, (37) 

from which follows that > w^. 

(ii) First, we will show that once a “stable system” has been reached, i.e. for fc > 0 we have 

P(w^) = the sequence will be constant. Let I be defined as above in (1321) . 

Note that such an I exists since and therefore P(w^"*“^ ) 5 P(w^) for all A: > 0 

as shown above. Due to r(w*) = r(w^+^), it follows because of 

^w*(w) = + M®r(w^)w = + M®r(w*^^)w = (w) (38) 

that the two mappings and are the same and consequently 

A direct consequence is that P(w*''“^) = P(w*'’“^) = P(w*) which implies r(w*'*“^) = 
r(w*"’'^) = r(w*). By iteration, all following vectors will be equal to 

What remains to be shown out is that the positive part of this iteration vector is the 
fixed point of the mapping <h®(-;f). Since is the fixed point of it holds that 

+ M®r(w*)w*+^. This yields to 

$s((w'+^)+; f) = (a + M'^f + M®(w'+^)+ - d)+ 

= (a + M'^f + M"r(w'+^)w'+^ - d)+ 

= (a + M'^f+ M"r(wV'^^-d)+ (39) 

= (w° + M®r(w')w''^^)’^ 

= (w'+')+. 

(iii) As shown, the series increases which means that the firms in P(w®) will maintain 

their positive entries in every further iteration step. The same statement holds for every 
firm i with < 0 and > 0 for any k > 0. Because of (ii) this means that the 

number of iteration steps would certainly be maximal, if in every iteration step the set 
P(w^) increased by one and if |P(w*"*“^)| = n. In that case we would therefore have 
|P(w^'’'^)| — |P(w‘^)| = n — do maximal possible iteration steps. 

□ 

Using Algorithm DAI to get an equity vector for a given debt payment vector, we can now 
present the algorithm to calculate the solution R*. In the sequel, we will make use of the 
mapping : (M))")"■ (M^)” defined by 

$‘^(r; s) = min{d, a + + M®s} (40) 

that represents the debt component of <I> for a given equity payment vector s > O^- 

ALGORITHM 3 (Elsinger Algorithm). Set e > 0. 

1. For k = 0, choose € {rsmaib^great} and determine s(r°) using Alaorithm \2A[ 
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2. For k >1, set ^;s(r^ ^)) and calculate s(r^) by Algorithmic 


3. If 


M—l 


< e, stop the algorithm. Else, set k = k + 1 and proceed with steplM 


The algorithm starts either assuming that all hrms can fully deliver on their debt obligations 
(rO _ = d) or that all firms have only their exogenous assets for paying their obligations 

(j,o _ = min{d, a}). With this payment vector, the corresponding equity payments are 

obtained by using Algorithm I2AI In the next step the debt vector has to be adapted to the 
new equity payments which is done applying to the previous debt vector. The updated debt 
payment vector is then used for determining a new equity payment vector. This procedure con¬ 
tinues until the i terates are suffi ciently close to each other. Additional to the original algorithm 
first presented in Elsinger ( 20091 1 . Algorithm [3] contains the second possible starting point rgmaii- 
We will show in the next proposition that if = rsmaii is chosen, the vector of debt and equity 
payments establish an increasing sequence and hence converges to the solution R* from below, 
while for = Tgi-eat; h converges from above. 

PROPOSITION 4. The Elsinger Algorithm delivers a series of decreasing vectors ifr^ = rg^eat 
and a series of increasing vectors if = rgmaii- Both series converge to the fixed point of the 
mapping <I> in dZD- 

Proof. The decreasing part is shown in Elsinger! ( 20091 1 . we only have to show that the debt 
iterate in the algorithm t herein is identic al to in Algorithm [HJ With our notation, the iterate 
of the debt component in Elsinger ( 20091 1 is defined as 


= min{d,(w*(r^-i) + d)+} 


(41) 


where w*(r^ is the solution of 


w = a + + M®w+ - d. 


(42) 


However, it follows from ([39]) that for r = 

(w'+i)+ = ((w'+i)+;r^-^) = (a + M^r^-i + M"(w'+^)+ -d)”^ = (w*(r^-^))^, (43) 

where is the result of Algorithm 12Al with the debt payment vector i.e. (w*(r^“^))'*' = 
s(r^“^). Because of (I42I1 . we have that 

w*(r^-i) = a + - d + M"(w*(r^-^))+, (44) 


from which follows with (14011 and a > 0,^ that 

=min{d,(w*(r^-^) +d)+} 

= min{d, a + + M"(w*(r^-i))+} 

= $d(rfc-i.(w*(rfc-i)) + ) 

= $‘^(r^-^s(r^-i)). 


(45) 


What remains to be shown is that for the starting point = rsmaii the generated series 
increases and converges to R*, which is done by induction. Eor the induction start check that 

= min{d, a} < min{d, a + -|- M®s(r°)} = 4>'^(r°; s(r*^)) = r^. (46) 
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As shown in Elsinger ( 20091 1 . the result w*(r) of Algorithm I2AI is increasing in r from which also 
follows that s(r) is increasing in r. Hence, s(r‘^) < s(r^) which completes the induction start. 
Assume for the induction step that and consequently s(r*^“^) < s(r^). The next debt 

iterate emerges as 

= min{d, a + + M^(s(r^))+} 


> min{d, a + + M"(s(r^-^))+} 


(47) 


= r 


from which also follows that s(r^'’“^) > s(r^) and, hence, the increasing property of the series. 
For the convergence, check that s(r^) > 0,^ and it holds that 

s(r^) = (a + + M®s(r^) - d)+ < M^s(r^) + (a + - d)+. (48) 

Because of < r*, it follows after some rearrangements that 

s(r^) < (I„ - - d)+ < (I„ - M^)-\a + - d)+, (49) 

hence the series s(r^) is bounded from above as well and therefore converges to some s* from 
below. The fact that 4>® is continuous in (r*,s*)* implies together with 4>®(s(r^); r^) = s(r^) 
that <l>®(s*;r*) = s*. Thus ((r*)*, (s*)*)* solves ([6]). Similarly, we can argue that because of 
the continuity of <h‘^(r*;s*) = r* from which follows that ((r*)*, (s*)*)* also solves dH) and 
therefore must be the fixed point R*. □ 


As described above, the Elsinger Algorithm determines the equity component of the iterates 
in a different way than the Picard Iteration. An important consequence of this approach 
is that the iterates of the Elsinger Algorithm will for the decreasing version be in every step 
smaller than the iterates of the Picard Algorithm, as we will show in the next proposition. 
The same statement holds for the increasing version of both procedures, where the iterates 
from the Elsinger Algorithm will be greater than the iterates form the Picard Algorithm. Both 
procedures are difficult to compare concerning their total calculation effort due to different ways 
of obtaining the next equity iterate (cf. Section (5.dp . However, if we only take the number of 
needed iterations as a quality criterion, we can conclude that the Elsinger Algorithm converges 
faster to R* than the Picard Iteration, no matter whether the algorithms start from the upper 
or the lower boundary. 

PROPOSITION 5. Let Rp = ((fp)*, (sp)*)* be the k-th iterate of the Picard Algorithm and 
R.| = ((r| )*, (sp)*)* the corresponding iterate of the Elsinger Algorithm. 

(i) For any iterate k > 0 it holds that Rp > Rp i/Rp = Rgreat Rp = (rgj,g 3 ^,-, s(rgreat)*)*- 
In case of = Rsmaii and R^ = (r Lai n s(rsmau)*)S we have that R| < R| for every 
iterate. 

(a) Let R^, k > 1, be an iterate either of the Picard Algorithm with R^ = Rgreat or of the 
Elsinger Algorithm with R° = (rgj,ga^^,s(rgreat)*)*- Then Rp'''^(R^) > Rp'^^(R^). If the 
starting vector is either liP = Rsmaii or R*^ = s(rarQaii)*)S it holds that Rp^^(R*^) < 

R^+^(R^). 

Proof. (i) The assertion is shown by induction. For k = 0, suppose that the upper boundary 
is the starting vector for both algorithms. In Equation (14911 it was shown that 

s| = s(d) < (In - M®)“^(a + M'^d - d)+ = Sgreat = Sp- (50) 
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Since = rp = d, the induction start is complete. Assume now, that for /c > 1 it holds 
that Rp > Rp. From Proposition 01 we know that Rp^^ < Rp. This leads to 

= min{d, a + + M®s^} > min{d, a + M‘^r| + M®4} = 4"^^ (51) 


= (a + M'^r^ + - d)+ 

> (a + + M"s| - d)+ 

> (a + - d)+ 

_ c^+l 

— »E ■ 


(52) 


If the starting vector is the lower boundary and the series Rp and Rp are increasing, the 
argumentation is similar. 


(ii) We prove the claim for the decreasing version of the algorithms, the proof for the reverse 
direction is similar. First, let R^ = Rp. The next iteration of the debt component is equal 
for both algorithms, i.e. rp"''^ = <h‘^(r^;s^) = rp"*"^. For the equity component, it holds 
that Sp"*"^ = <I>®(s^; r^). The mapping r*') has a unique fixed point, that we denote by 
s(r^) and that can be obtained via a Picard Iteration: 


lim($>)' (s‘;r'')=s(r'=). 

/—>-oo 


(53) 


The iterates obviously form a decreasing sequence so that 






> s(r'=) > s(r"+^) = (54) 

where the second inequality follows from the fact that s(r) is increasing in r (cf. Elsingeil 
(2009|)). If the k-th. iterate is given by R*' = Rp, the arguments are analogous to the ones 
above. 

□ 


3.3. A Hybrid Algorithm 


To motivate the approach of the next algorithm, we have to compare the functioning of the 
Elsinger Algorithm and the Picard Iteration. The major difference between both iterations 
emerges in the calculation of the equity component. Suppose that we are in iteration step 
A: > 0 and want to calculate the next iteration of the equity component. We ignore for an 
instant that both algorithms deliver different iterates and assume that the fc-th iterate is given 
by R*^ = ((r^)*, (s^)*)*. In the Elsinger Algorithm, is calculated first and then as the 
fixed point of so that it holds that r^+^). The Picard iterate, on 

the other side, can be written as = 4>®(s^;r^) from which it becomes clear that the Picard 
Iteration neither uses the “updated” debt vector nor does it solve a separate fixed point 

mapping to obtain 


The determination of the debt component 




however, is comparable in both algorithms. 


Again starting with R^ we have that = <I>‘^(r^;s" 


for both procedures. An obvious exten¬ 
sion of the Elsinger Algorithm would be to utili ze the principle used for t he equity component 


for the debt component as well. In the article of Eisenberg and ( 20011 1. this concept is used 
for systems with no cross-ownership of equity, i.e. where M® = 0„xn- In this subsection we will 
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generalize the results of this work and it will turn out t hat combining both ideas, the one of 
Elsinger ( 20091 ) and the one of Eisenberg and Noe (2001), will help to minimize the number of 


needed iteration steps of the global algorithm to find R*. 

To explain this idea in more detail, say that for a debt payment vector f G [rsmall j Tgreat] we 
have a corresponding equity vector s(f), that is, a fixed point of the mapping <l>®(-;f) in (1261) . 
In the Elsinger Algorithm the next debt iterate emerges as s(f)). Instead of using this 

iterate, our aim is now to find the fixed point of s(f)) as the new iterate. This can be done 
using the following Algorithm. 


ALGORITHM 4A. Suppose s > 0„. 

1. For k = 0, set = f and determine D{r^,s) and A(r^,s). 

2. For k > 1, solve ^ where 

Q,.k-i s{r) = A(r*^“^,s) (^A(r^ 

+ - A(r*^"\s)^ d 


-1 


r+ (R-A(r^-\s)) d) +M®s 


(55) 


3. Denote the solution by , i.e. 0j.fc-i s(''^) = determine D{r^,s) and A(r^,s). 

4- If D{r^,s) = ZI(r^“^,s), stop the algorithm. Else, set k = k + l and proceed with step\^ 


The algorithm is identical to the one given in lEisenberg and Noel (j200ll ) with the modification 
that some additional fixed payments due to equity cross-ownership are included. It solves ([5]) 
for a fixed amount of equity payment s > 0„, i.e. is the fixed point of the mapping <l>‘^(-,s), 
as we will show in the next proposition. Denote this fixed point by r*(s) for instance. In the 
Hybrid Algorithm following later on, r*(s) is used as the next iterate for the debt component. 
To see the difference between the calculation of the debt component in the Elsinger Algorithm, 
assume that an arbitrary debt payment vector r G [0„,d] is given and that the corresponding 
equity payment vector s(r) is given too. The fixed point of the mapping <l>^(-;s(r)) can on the 
one hand be obtained using Algorithm I4AI above, but on the other hand, we could also use a 
Picard Iteration, since for any r G [On,d] it holds that 


On < s(d)) = min{d, a -|- M^r + M®s(d)} < d. 


(56) 


Starting with the vector r, the fixed point r*(s) is given as 

r*(s) = lim (<l>‘^)*(r; s(d)). (57) 

l—^OO 

In the Elsinger Algorithm, however, the next iterate for the debt component is defined as 
‘h‘^(r; s(r)) which is therefore the first iterate of the Picard Iteration in (I57p . Hence, one can say 
that using in the Elsinger Algorithm, a simple mapping is applied to obtain the next iterate, 
whereas in the Hybrid Algorithm, the fixed point of a mapping is determined. In Proposition 
[9l we will show that when using the idea of the latter algorithm, the iterates R^ will always be 
closer to the searched solution R*. 


PROPOSITION 6. Let f G 

equity payments such that 


[rsmaii) Tgreat] ^6 a debt payment vector and s > 0„ 

‘h‘^(r; s) < f. 


a vector of 


(58) 
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(i) Algorithm \4A\ generates a well-defined decreasing sequenee of vectors r^. 

(ii) Let 1 <l <n such that 

I := min{j € {0,1,... ,n} : D{rfis) = s)}. (59) 

Then r*(s) = is the fixed point of the mapping defined in iO]). 


(Hi) Let do = |D(r, s)| be the number of firms in default under r and s. The fixed point r*(s) 
is reached after no more than n — do iteration steps. 

Proof. Since the equity vector s is considered as fixed we can modify the financial system P by 
setting a = a + M®s and M® = Onxn- The new system T = (a, d, M®) is then a syst em 
without cross-ownership of equity. Such systems are considered in Eisenberg and ( 200 li b 

(i) The proof that the seq uence decreases is now equivalent to the proof given in 
Eisenberg and ( 2001 b A needed assumption in the proof therein is that f is a so- 


called supersolution which is given because of <l>‘^(f; s(f)) < f < d. What we have to 
show to complete this part is that the fixed point of the mapping in (l55]( exists and is 
unique, since their definition of a financial system, differs slightly from ours. Denote by 
A := A(r^,s) the diagonal matrix for r^. The next iterate is according to (I55p given 
by 

A+i = A (a + M'^(Ar^+^ 

(60) 


+ {In - A)d)) + {In - A)d 

= AM‘^Ar^+^ + A (^a + M‘^(I„ - A)d) + (I„ - A)d 
and rearranging yields to 


= 


In - AM^A'j ' (a (a + M^{ln - A)d) + (I„ - A)d) . 


(61) 


Note that has the Elsinger Property and, hence, so does AM'^A, which means that 
the inverse of I„ — AM'^A exists. This proves the uniqueness of 


(ii) The argumentation that the sequence converges and becomes constant in the end is 
analogous to part (ii) of the proof of Proposition [3l Since is decreasing, we have 

that D{r^,s) C T)(r^'*“^,s) that means the number of firms in default increases. If 
D(r^,s) = D(r^+^,s), then we also have that A(r^,s) = A(r^+^,s), from which follows 
that the mappings 0j.i g and 0j.i+i g have the same fixed point. It must hold then that all 
consequent iterates are equal. 

To show that is the fixed point of ‘h‘^(-,s), first check that by definition of A(r^'^^,s): 

r'+i = A(r'+\ s)r'+^ -h (In - A(r'+\ s))d. (62) 

It then holds that 

$‘^(r'+T s) = min{d, a + 

= (In - A(r'+\ s))d + A(r'+\ s)(a + M^r'+i) 

= (I„ - A(r',s))d 

-h A(r',s) (^A(r'+\ s)r'+^ -h (In - A(r'+\s))d^^ 

= (In - A(r',s))d A(r',s) (^a(^A(r', s)r'+^ -h (In - A(r',s))d^^ 

= 

(63) 
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where the last equality follows from (155 p . 

(hi) This part is similar to part (hi) of the proof of Proposition [3] with the reverse argumen¬ 
tation. The do firms in default under the starting vector will stay in default since the 
series decreases. To achieve a maximum theoretical length of the algorithm, exactly one 
additional default step has to occur in every new iteration step. This results in no more 
than n — do possible iteration steps. 


□ 

The validity of the inequality in (|58p is crucial for the monotonicity of the iterates produced 
by Algorithm I4A1 However, there are situations in which a debt payment vector f € [rsmaib ^great] 
is given together with an arbitrary vector s > 0,^ and where (I58p does not hold. Think of an 
algorithm to find R* that starts with Rsmaii- In this case the first debt iterate is rgmaii = 
min{d,a} and the corresponding equity iterate is s(rsmaii)- Applying on these vectors yields 
to 


^'‘^(rsmali; s(rsmall)) = min{d, a -f M'^rs^all + M®s(rsmall)} > min{d, a} = rgmall (64) 

and to a violation of (|35|1 . Finding the next debt iterate as the fixed point of s(rsinaii)) 
and applying Algorithm I4AI to do so, can under certain circumstances lead to a non-monotone 
series, as one can simply verify by a self-chosen example. This makes it difficult to prove the 
convergence of such a series in general. Nevertheless, given a debt vector f and s > 0,^, we can 
still calculate the fixed point by avoiding Algorithm 14Al and use a Picard-type algorithm instead. 

ALGORITHM 5A (Picard Iteration for the Debt Component). Suppose that s > 0„ and 
e > 0. 


1. For /c = 0, set = r. 


2. For k>l, determine ^;s). 

3. If ||r^“^ — r^ll < e, stop the algorithm. Else, set k = k + 1 and proceed with stepO 

PROPOSITION 7. Algorithm\M delivers a series of decreasing vectors if <I>‘^(r;s) < r 
and a series of increasing vectors i/<h‘^(f;s) > f. Both series converge to the unique fixed point 
of^%;s). 


Proof. Assume that <I>‘^(r; s) > r = r°. For the first iterate, it holds that = <I>'^(r; s) > Via 


induction, it follows that for all k > 1. 


Because the monotone series is bounded 


by d, it must converge to some fixed point. Because of th e fact that the Elsinger condition 
holds, it follows directly that this fixed must be unique (see Elsinger ( 20091 1. Theorem 3). The 
argumentation is similar if <h‘^(f;s) < f. □ 


The Algorithms I4AI and DAI both enable us to calculate a new debt iterate given an equity 
vector. Together with Algorithm DAI for the equity component, we can now combine both 
procedures in a common algorithm that searches for the fixed point R*. 


ALGORITHM 6 (Hybrid Algorithm). Set e > 0. 

1. For k = 0, choose G {rgreat,rsmall} CLnd determine s(r‘^) with Alaorithm \2A[ 

2. For k > 1: 
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2.1 Determine using Algorithm \4.A\ if = fgreat or using Alaorithm, \5A\ if = rs mall 
in both cases with s = s(r^“^). 

2.2 Determine = s(r^) using Alaorithm \2A[ 

\ 

I < e, stop the algorithm. Else, set k = k + 1 and proceed with step\^ 

For given r = r^,A: > 0, the Hybrid Algorithm determines = s(r*^) as the correct equity 
value that solves ([6]) and for given s = s(r^), A: > 0, it determines the correct debt value that 
solves ©. As such, conditional on the values determined in the previous step, the algorithm 
calculates an exact solution of either ([6]) or ([5]) in the next iteration step. 

PROPOSITION 8. The Hybrid Algorithm delivers a series of decreasing vectors i/r° = rgreat 
that converges to the fixed point R* . In case of = rgmaii the series is increasing with the same 
limit. 



Proof. First, suppose that = Ugreat- We will first show by induction that the series decreases. 
For the induction start note that 


s(r*^)) = min{d, a + + M®s(r°)} < d = r°. (65) 

As mentioned in the proof of Proposition O the equity vectors s(r) are increasing in r which 
yields to s(r^) < s(r‘^). For the induction step, assume that for /c > 1 it holds that 
and consequently s(r^“^) > s(r^). Since s(r^“^)) and because of 

^ min{d, a + + M®s(r^)} 

< min{d, a + + M"s(r*^-^)} 


the assumption (|58p is fulfilled. The next iterate emerges from a decreasing sequence 
produced by applying Algorithm l4Al beginning with f = r^. Hence j,fc+i < j.fc < 

s(r^). Next step is to show that the series converges to R*. We have that the two sequences 
((rA:+i)t, (s(r^))*)* and ((r^)*, (s(r^))*)* are both decreasing in and therefore converge to 

the same limit ((r*)*, (s*)*)* € Because of the continuity of and it must hold 

that <l>‘^(r*,s*) = r* and <l>®(s*,r*) = s*. Thus, ((r*)*, (s*)*)* solves ([6]) and ([5]). The proof for 

= Tsmall is similar. □ 

In Proposition [5l we have shown that when using the Elsinger Algorithm, the iterates will 
always be nearer to the solution R* than the corresponding iterates of the Picard Algorithm. 
This lead to the conclusion that the iteration number is minimized for the Elsinger Algorithm. 
The next Proposition shows the same when comparing the Elsinger and the Hybrid Algorithm 
and it will become clear that the Hybrid Algorithm will need less iteration steps to reach R* 
than the Elsinger Algorithm. 

PROPOSITION 9. As in Proposition O we denote the iterates of the two algorithms with 
subscripts, where E stands for the Elsinger and H for the Hybrid Algorithm. 

(i) For any iterate k > 1 it holds that Rg > R^ if = (rgj.g 3 ^^, s(rgreat)*)* o.nd R^ < R^ 
when (rgrrjg^ii, s(rgjnau)*)* is the starting vector of both algorithms. 
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(ii) Let R^, A: > 0, be an iterate either of the Elsinger Algorithm or of the Hybrid Algorithm 
that started with = (rgj.g 3 ^^, s(rgreat)0*- Then Rg^^(r^) > R^’^^(r^) for the next iterates 
which were calculated with either the Elsinger or the Hybrid Algorithm starting from R*^. 
//R° = ((r Lai bS(rsmaii)*)S it holds that < R|+^(r''). 

Proof. (i) Let R° = (rgj.g 3 ^^, s(rgi.eat)*)*- From Proposition [8] we know that = d which 

yields to 

= min{d, a + M^d + M®s(d)} > min{d, a + M'^r^ + M®s(d)} = r^. (67) 

Further, since s(r) is increasing in r (cf. Proposition [5|), = s(rg) > s(r^) = s^, which 

completes the induction start. For the induction step, assume that it holds for /c > 1 that 
Because of Proposition [HI and thus 

r| = min{d, a + + M®s(r^-^)} 

> min{d, a + + M®s(r|-^)} 


R k—l \ ■D ^—1 

E — -^H • 


> min{d, a + + M"s(r^-^)} 


„k-l\ 


( 68 ) 


— 




where we again used the fact that s(r) is increasing in r from which follows that s(rg j ^ 


> 


s(r^ ) and also s(rg) > s(ry). The proof when R° = (r smai u^(^^smaii)*)^ is completely 


k-l 


analogous. 

(ii) Let R*^ = (rgj,ga^,-, s(rgreat)*)* and R^ = R^. Note that because of 

^>‘^(r|; s(r|)) = min{d, a + + M^s(r|)} 

< min{d, a + + M®s(rg“^)} 

- 

— ^E- 

Therefore, the assumption in (I58p is fulfilled which ensures that < r 
iterate it follows that 


is holds that 


(69) 


= min{d, a + + M^s(r^)} 

> min{d, a + + M"s(r|)} 


For the next 


(70) 


_ ^fc+i 


= r 


which in turn implies 

H ^ ^ H to 

/c+1 _ 


H 

/c+1 


because of < r^j to 


On the other hand, starting with R^ = R^ yields 


= min{d, a + + M"s(r^)} 


> min{d, a + + M®s(r^)} 


(71) 


= r 


fe+i 

H • 


If follows from this results that A similar argumentation together with 

Proposition 13 delivers the proof in case of R° = s(rsniaii)0*- 

□ 


Of course, within an iteration step of the Hybrid Algorithm, potentially many linear equation 
systems have to be solved since for the debt component Algorithm I4AI is applied, which results 
in higher computational costs. But if we ignore for a moment this circumstance it follows from 
Proposition [9] that the convergence speed of the Hybrid Algorithm is higher than the one of the 
Elsinger Algorithm. 
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4. Finite Algorithms 


The Algorithms in the previous section all had the drawback that it could not be ensured that 
the solution R* is reached in a finite number of iteration steps. In this section we will present 
two ways in which potentially infinite solution algorithms can be turned into procedures that 
reach the solution in hnitely many steps. The common principle of these methods is to include 
the information which firms are in default under a current iterate It turns out that this 
slight modification helps to overcome the disadvantage of potentially infinitely many iteration 
steps. 

To guarantee that the forthcoming procedures are well-defined, we ha ve to drop the Elsinger 
Property and demand a stricter property of the ownership matrices (see [FischeiJ (|2ni4l ')'l. 


ASSUMPTION 2. For both the debt and the equity ownership matrices it holds that HM'^H < 1 
and ||M®|| < 1. 


For the remainder of this section we suppose that Assumption [2] holds. Note that Assumption 
[2] implies Assumption [H but not the other way round. The financial system therefore still has 
a unique solution under Assumption [2l 


DEFINITION 4. Let R = (r‘ 

set D(R) and default matrix A 
belongs to D(R) is defined by 


, s*)* € be an arbitrary vector with corresponding default 

= A(R). The pseudo solution R € of (l5|) and ([6]) that 


AU-A)d + Ax\ 

V (In-A)x )' 


(72) 


where x G M” is the solution of the linear equation system Ax = b with 


A = ln- (m'^A + M"(U - A)) G 


(73) 


and 

b = a + (U - A)d - (U - A)d G R”. (74) 


To motivate the definition of a pseudo solution, assume that it was known for each firm 
whether it was in default under the solution of ([5]) and Q or not. Denote by D* C J\f the set 
of firms that were in default under R*: 


+ p^M!jS*<d)l (75) 

and let A* = A(r*, s*) be the corresponding default matrix. We assume that the set was known 
even though this information is not available a priori. However, if we had this information, no 
iteration procedure would be needed to find the fixed point R*. We only had to compute the 
pseudo solution that belongs to D*, as is shown in Proposition [TUI 

The reason why we have to restrict the following considerations to ownership matrices with a 
matrix norm smaller one is because we have to guarantee that x from Definition [4] is uniquely 
defined. This can only be ensured if II < 1 and ||M®|| < 1 since then ||M^A-|-M®(U —A)|| < 
1 for any A as well. This in turn implies that A in (17311 is invertible. For ownership matrices 
and M® that have the Flsinger Property, the invertibility of A is obviously not always given. 


f " 

D* = D{r*,s*) = ii€M:ai + Yl 
[ i=i 
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PROPOSITION 10. The pseudo solution belonging to D* is the solution R* of the financial 
system iF {a., i.e. 


R* = 


(R - A*)d + A*x' 
(R - A*)x 


(76) 


where A* is the default matrix belonging to D* and x is the solution of the equation Ax = b 
defined in (I73]l and 

Proof. According to the liquidation value equations in (l5|) and ([6|), the vectors r* and s* are 
given as 

^ fdj, if i ^ D*, 


(77) 


and 




0 , 


if i G D* 


(78) 


R* = 


(79) 


In matrix notation this means in particular that (I„ — A*)r* = (I„ — A*)d and A*s* = 0^ and 
thus 

'(R-A*)d + A*r* 

(In-A*)s* 

For the firms in default we only have to calculate the debt payments and for the firms not in 
default we have to determine the equity value. The solution R* does hence contain only n 
unknown values and we only have to consider the two subsystems 


A*r* = A*a + A^M'^r* + A*M®s* 


and 


(R - A*)s* = (R - A*)(a + + M®s* - d) 

We can add the two equations and write the system more compact as: 

A*r* + (R - A*)s* = a + + M®s* - (R - A*)d. 

Because of (!„ — A*)s* = s* we get 

A*r* + (R - A*)s* = a + + M^(R - A*)s* - (R - A*)d, 

which leads after some rearrangements to 

A*r* + (R - A*)s* - M'^AW* - M®(R - A*)s* = a + M‘^(R - A*)r* - (R - A*)d (84) 


(80) 

(81) 

(82) 

(83) 


that is equivalent to 

(r- (m‘1A* + M^(R-A*))) (AW* + (R-A*)s*) = a + M‘i(R-A*)d-(R-A*)d, (85) 

since A*(I„ — A*) = 0„xn- Setting x = A*r* + (I„ — A*)s* and with the notation of Definition 
m the equation system becomes Ax = b. □ 


The main challenge in this solution approach is of course that the final default set D* is 
unknown. Algorithms that follow this idea to find R* consequently have to find D* in a fast way. 
A naive strategy could be to check all possible default scenarios of the financial system, calculate 
the pseudo solution for the corresponding default set and check whether it actually is the fixed 
point of However, there are 2"’ possible scenarios that would have to be checked, which could 
be cumbersome for large n. Therefore, more efficient algorithms are needed that require less 
computation to find D*. Some possible algorithms are presented in the next subsections . 
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4.1. Decreasing Trial-and-Error Algorithms 

The three algorithms [U [3] and E] from Section [3] can start with a vector that was the upper 
boundary of the solution vector R*. The procedures in this subsection have in common that 
they also start with this upper boundary and calculate a corresponding default set. For every 
following iterate, the corresponding default set is determined as well. To avoid that every default 
set it is checked whether it actually is D* and whether the corresponding pseudo solution is the 
fixed point of the algorithm will identify potential default sets to reduce the computational 
effort. If it turns out that the potential default set is D*, the algorithm stops. Otherwise, the 
procedure continues until a new potential default set is found that has to be checked again, and 
so on. Due to these characteristics we name this type of algorithm Trial-and-Error Algorithm. 
The general procedure of algorithms of this type is similar. 

ALGORITHM 7 (Decreasing Trial-and-Error Algorithm). Set I > 2 and p = 0. 

1. Choose either the Picard (Algorithm\J^, the Elsinger (Algorithmic or the Hybrid Algorithm 
(Algorithmic which is used in the following to generate the next iterate. 

2. If in Step [7] the Picard Algorithm is chosen, set d = —1, R^ = Rgreat an-d determine 
T)(R‘^). Else, set d = 0, R*^ = ^s{d) ) determine DfRP). 

3. //D(R0) =M, setK* = and stop the algorithm. 

4 . If the Elsinger or the Hybrid Algorithm is chosen in Stepl^and if = 0, set R* = R*^ 

and stop the algorithm. 

5. Else, calculate for k > p the iterates R^ starting with R^ using the algorithm chosen in 
step [7] and the corresponding default sets D(R^) until k = q with 

q = min{m > p : = ... = D(R'") and |D(R™')| > d} (86) 

is reached. Determine the pseudo solution belonging to D(R'^) and denote it by R*?. 

6. //<1>(R'^) = R”?, stop the algorithm. Else, set d = |D(R'^)| and p = q and proceed with step 

m 

The Algorithms [Ills] and [6] in their decreasing versions produce decreasing sequences of iterates 
and thus increasing sequences of default sets, i.e. D(R^) C D(R^"''^) for A: > 0. Algorithm [7] 
means that one iterates and checks whether the default set has not changed compared to the 
previous default set. If the default set stays the same for the next I — 1 consecutive iterations, 
this is an indication that the actual D* might have been reached. To check this, the pseudo 
solution is calculated and it is checked whether it solves © and (E|). If no solution has been 
found, one iterates again until a larger default sets stays identical for I — 1 consecutive times, 
and the described procedure can be repeated. If a solution is reached, the procedure stops. Due 
to its described property, we call I the lag value. 

In the special case of I = 2 this means that the pseudo solution is calculated if the default set 
stays the same from one iteration step to another. Obviously, choosing a higher lag value inspires 
more confidence in the potential default set since the longer the default set stays unchanged, 
the higher is the chance that it is the actual default set. 

Depending on the choice of the algorithm in Step [T] of the Decreasing Trial-and-Error Algo¬ 
rithm, we obtain three different versions of Algorithm [7) 
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(i) The Decreasing Trial-and-Error Picard Algorithm with = Rgreatj where the iterates as 
given by R^ = <h(R^“^). 

(ii) The Decreasing Trial-and-Error Elsinger Algorithm with R^ = ((rgreat)*, (s(rgi-eat))*)*) 
where s(rgreat) is obtained via Algorithm I2AI and the next iterates are obtained using 
Algorithm [3l 

(iii) The Decreasing Trial-and-Error Hybrid Algorithm with the same starting vector as in (ii), 
where the next iterates are obtained using Algorithm [6l 

The particular cases when T)(R‘^) € {0,A/'} in the steps [3] and U deserve a separate mention 
since in such situations, no iteration is necessary and the solution R* can be given explicitly under 
some circumstances. The justification of this phenomena is given in the following proposition. 

PROPOSITION 11. Eor the Decreasing Trial-and-Error Hybrid Algorithm the following 
holds: 

(i) //T>(R°) = M, then R* = ( matter which version of the algorithm is 
taken. 

(ii) If Z1(R0) = 0 and either the Decreasing Trial-and-Error Elsinger Algorithm or the De¬ 
creasing Trial-and-Error Hybrid Algorithm is used, then R*^ = R*. 

Proof. (i) First, assume that the Picard Algorithm is chosen in Step [T] of Algorithm [71 Because 
of DfRP) = J\f, it must hold that a + M'^d + M^Sg^eat < d and also a + M'^d < d. A 
consequence is that Sg^eat = (In — M®)“^(a + M'^d — d)+ = 0„. From Proposition [1] it 
follows that s* = 0. For s = s* = 0, Equation ([5]) is now solved by r* = (I„, — 
where Lemma [A3] proves that (I„, — exists. If the Elsinger or the Hybrid Algorithm 

is chosen in Step [U we have that a + M'^d + M®s(d) < d. It follows from ([77|) that 
s(d) = s* = 0,1 since s* = s(r*) < s(d) because of r* < d and the fact that s(r) is 
increasing in r. The solution of Equation ([2]) is therefore the same as in the Picard case. 

(ii) Now, R° = (^s(d)) since T>(R°) = 0, it holds that a + M‘^d + M®s(d) > d. This leads 
to 

(hf d A /min{d,a + M‘id + M"s(d)}\ / d\ 

Vs(d)y V(a + M‘id + M®s(d)-d)+y vs(d)y’ ^ 

which proves the claim. 

□ 

Note that for the Decreasing Trial-and-Error Picard Algorithm, we cannot conclude that 
R° = Rgreat = R-* if Il(R^) = 0. There are simple counterexamples for situations like this. 

PROPOSITION 12. Algorithm^ reaches the solution R* of ([5]) and ^ in a finite number 
of iteration steps. 

Proof. By definition of DfR^) in (1141) and since R*' converges to R* from above for any of the 
three algorithms [H [3] and [H there exists a > 0 such that Zl(R^) = D(R*) = D* for all 
k>k^. □ 


24 


4.2. Increasing Trial-and-Error Algorithms 

In contrast to the decreasing algorithms presented in the subsection above, it is of course also 
possible to use an algorithm with the reverse direction, i.e. in which the series of produced 
iterates is increasing and in which the default sets are decreasing. The general form is very 
similar to Algorithm [71 

ALGORITHM 8 (Increasing Trial-and-Error Algorithm). Set l>2, d = n + l and p = 0. 

1 . Choose a starting vector R*^ and determine H(R^). 

2. If D(RP) = 0, set R* = ^s(d)) the algorithm. 

3. Else, calculate for k > p the iterates R^ starting with R^ using one of the Algorithms [21 
0 or0 and the corresponding default sets DfR^) until k = q with 

q = min{m > p : T)(R™“'+^) = ... = DfR^) and |T>(R™)| < d} (88) 

is reached. Determine the pseudo solution belonging to Z?(R'?) and denote it by R*?. 

4. //4>(R5) = R'', stop the algorithm. Else, set d = |H(R'^)| and p = q and proceed with step 

0 

The functioning of Algorithm 0 is similar to the Decreasing Trial-and-Error Algorithms with 
the difference that the resulting sequence of default sets is obviously decreasing. As in Section 
HU the way of choosing the calculation method to determine the next iterate, allows three 
different modifications; 

(i) The Increasing Trial-and-Error Picard Algorithm with R^ = Rgmaii and R^ = <I>(R^“^). 

(ii) The Increasing Trial-and-Error Elsinger Algorithm with R*^ = ((rgmaii)*) (s(rsmaii))*)*) 
where s(rsinaii) is obtained via Algorithm I2AI and the next iterates are obtained using 
Algorithm 0 

(iii) The Increasing Trial-and-Error Hybrid Algorithm with the same starting vector as in (ii) 
and where the next iterates are obtained using Algorithm 0 Note that for the next debt 
iterate, Algorithm I5AI is used instead of Algorithm I4AI in the decreasing version. 

The justification of the stopping criteria in Step [2l of Algorithm 0 is as follows. Suppose 
that the Picard version of the algorithm is chosen and that Z)(R^) = 0, which means that 
a -|- IM Tgujaii T NI Ssmaii ^ d. Since rgmaii < d and Sgmaii < s(d), it also holds that a -|- -|- 

M®s(d) > d and s(d) > 0^ following from this. With Equation (l87|) . we see that R* = (^s(d)y 

Also note that, in contrast to Algorithm 0 there is no stopping criteria in case of D(R^) = J\f. 
The reason is that in this case, no general statement about the structure of the solution R* can 
be made, no matter which version of the algorithm is used. In particular, from ZI(R^) = A/" it 
does not follow in general that D(R*) = Af, since there are easily constructable counterexamples 
for this. 

The reason why we distinguish between decreasing and increasing Trial-and-Error Algorithms 
is that Algorithm 0 will always find the correct default set D* = D(R*), and this in a finite 
number of iteration steps. For the Increasing Trial-and-Error Algorithms such a statement is not 
possible in general since there are some situations in which the default sets do not converge to 
D*, no matter which lag value is chosen. Situations in which this “anomaly” occurs are always 
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finan cial systems that cont ain a so-called borderline firm. The expression borderline is taken 
from Liu and Staum ( 2010l l and denotes a firm i G AA in a financial system with fixed point R* 
for which it holds that r* = di and s* = 0. In other words, borderline firms are just able to 
fully cover their liabilities, but have no remaining capital left in their balance sheet that can be 
furnished to their shareholders. By definition of a default set in (|14D . a borderline firm i is not 
in default since 


Q = Si = ai + Y,Mtf* + Y.M\ 
i=i i=i 


* /7 


(89) 


and therefore i ^ D(R*). However, when using an Increasing Trial-and-Error Algorithm it can 
happen for such a borderline firm i that i € I?(R^) for every iterate R^, k > 0. This means that 
the true default set D* will never be identihed by the algorithm. There exist many examples of 
financial systems that have this property. To show that in such situations, the fixed point R* 
can still be determined via the calculation of the pseudo solution, assume that the set B C Af 
contains an arbitrary selection of borderline firms. The common set of defaulting firms and 
the selected borderline firms is denoted by D, i.e. D = D* U B. The corresponding “default” 
matrices are given by A = A.[D) and A* = A*{D*), respectively. Following this notation, A 
and A* define the matrices from ()73p with the corresponding default matrix, and b and b* are 
defined analogously. Moreover, we define Ag = A — A* as the diagonal matrix that indicates 
only the selected borderline firms. 


LEMMA 2. The vector x* = A*r* -|- (R — A*)s* solves the equation system A*x = b* if and 
only ifx = x*+ Agd is the solution of Ax = b. 


Proof. Without loss of generality, we assume that the hrst ni hrms of the system are solvent, 
that the next n 2 — ni firms are the selected borderline cases, and that the remaining firms are 
in default under R*. This means that 


Af ={!,...,ni}uBuD* = {!,.. 

It follows from Proposition 1101 that 

., ni} U {ni + 1,..., 77 - 2 } U {n2 + 1,... 

,n}. (90) 

X — {xi, ... , Xni , Xni+1 j ■ ■ ■ j Xn2 j Xn2+1 j ■ 

Further, note that 

..,XnY = (st,...,s*^,0,...,0,<2+i,. 

..,<)*. (91) 

M^(I„ - A*)x* = M^(In 

and 

-A)x* = M^(R-A)(x* + Aed) 

(92) 

M'^(I„ - A*)d - 

and that 

M'^Agd = M^(R - A)d 

(93) 

MdA*x* + M^Asd = M‘^A(x* + A^d) 

(94) 


because of the structure of x*. By d73|) and (ITill . x* solves A*x = b* if and only if 


X* = b* + (M^A* + M^(I„ - A*))x* 

= a + M'^(I„ - A*)d - {In - A*)d + (M^A* + M®(I„ - A*))x* 

= a + M'^(I„ - A*)d - M'^Agd - (I„ - A*)d + (M'^A* + M"(R - A*))x* + M'^Afjd 
= a + M^(R - A)d - (In - A*)d + (M'^A + M^(R - A))(x* + Agd). 

(95) 
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Since (I„ — A*)d = (I„ — A)d + Agd, we can add Agd on both sides of the equation and obtain 

X* + Afjd = a + - A)d - (I„ - A)d + (M'^A + M^(I„ - A))(x* + Agd), (96) 

Therefore x = (s*,..., s*^, dm+i, ■ ■ ■, dn 2 ,r*^+i, ■ ■ ■, is the solution of Ax = b if and only 
if X* solves A*x = b*. □ 

The pseudo solution belonging to D* is the solution R* of the system. A direct consequence of 
Lemma[2]is that the pseudo solution of D is also equal to R*. Similar to the proof of Proposition 
m we can argue that the set D will be reached by the Increasing Trial-and-Error Algorithms in 
a finite number of steps. Note that this statement holds in particular for the Increasing Hybrid 
Trial-and-Error Algorithm, where Algorithm is used 15 Al to calculate the next debt iterate. Even 
though a Picard-typed procedure is used in this auxiliary algorithm, we can conclude together 
with Proposition [9] and e > 0 that the number of iteration will still be finite. We summarize the 
findings in the next proposition. 

PROPOSITION 13. Algorithmic reaches the solution R* of ([5|) and ^ in a finite number 
of iteration steps. 

4.3. Sandwich Algorithms 

A disadvantage of the Trial-and-Error Algorithms was that when a potential default set is 
reached, the only way to find out whether this default set is actually D*, is to calculate the 
corresponding pseudo solution and check whether it is a fixed point of ([7]). The choice of a high 
lag value can of course increase the chance that D* is reached at the first trial, but there is no 
certainty. 

Another way to find D* is to start an iteration simultaneously with the largest and smallest 

possible solution and use one of the Algorithms [H [3] or [6] to obtain the next iterate. For k > 0 

_ ^ 

denote by R the fe-th iterate of the series that emerges when starting the algorithm with 
the maximum and by R^ its counterpart when starting with the minimum possible solution. 
Depending which algorithm is chosen, the starting vector can either be R = Rgreat (Picard 
Iteration) or R*^ = (rgj,g 3 ^^-, (s(rgi.eat))*)* (Elsinger and Hybrid Algorithm) when starting with 
the upper boundary. Analogously, we have R*^ = Rsmaii or R*^ = (s(rsinaii))*)* if the 

minimum possible solution is the starting point. By the Propositions [21 01 [8] and by Equation 
m, the iterative use of one of the mentioned algorithms entails that the default sets approach 
one another, i.e. for k >0 

D{RY) D D(R^+1) T D* T T)(R^^^) T D(R^). (97) 

Let 

I = mm{k > 0 : D(R^) = D(R^)} (98) 

be the first iteration step in which the default set for both starting vectors is the same. Then we 
must have that DfR ) = D* and, by Proposition 1101 determining the pseudo solution belonging 
to D* leads to R*. Because of its characteristics we call this algorithm the Sandwich Algorithm. 

ALGORITHM 9 (Sandwich Algorithm). 

1. Determine ^ and R*^ as well as their corresponding default sets D(^) and DfRf). 
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_ ^ , 

2. For k>l, calculate the iterates R and R* using one of the Algorithms{^ 0 or0 and the 

_ ^ , 

corresponding default sets DfR ) and DfRf). 

3. If DfR^) = DfR!^), stop the algorithm, set D* = D(R^) and calculate the pseudo solution 
that belongs to D* following Definition^ Else, set k = k + 1 and go back to steplM 

As for the Trial-and-Error Algorithms in the sections above, the Sandwich Algorithm results 
in different versions: 

(i) The Sandwich Picard Algorithm with R*^ = Rgreat and R° = Rsmaii and the use of Algo¬ 
rithm [U in step [2j 

(ii) The Sandwich Elsinger Algorithm with R° = (s(rgreat))*)* and R*^ = 

(^smaiu (s(rsmaii))*)* and the use of Algorithm [3] in step El 

(iii) The Sandwich Hybrid Algorithm with the same starting points as the Sandwich Elsinger 
Algorithm and the iterative use of Algorithm [6] in step El 

Recall the insights from Section 14.21 where it was shown that, under some circumstances, 
it may happen that the series of default sets Tl(R^) will never converge to the actual default 
set D*. Situations in which this problem occurs always contain at least one firm that is on 
borderline in the solution R*. As a result of this behavior, the Sandwich Algorithm may not 
converge in the sense that the default sets D(R ) and D(Rf) will never become identical. 
However, if we consider a stochastic setting and assume a distribution for the vector a of the 
exogenous assets’ prices which has a density with respect to the Lebesgue measure on 
then situations in which the convergence cannot be assured occur only with probability zero as 
the next Proposition shows. Note that this assumption is fulhlled in the usual n-hrm Merton 
models where the individual a* are log-normally distributed. 

PROPOSITION 14. The Sandwich Algorithm generates a sequence of decreasing default sets 
T>(R'‘) and a sequence of increasing default sets Z)(R ) that reach the default set D* of the 
solution R* almost surely after finitely many steps. Thus, it reaches the solution R* of d?]) 
almost surely after finitely many steps. 


Proof. The increasing and decreasing property of the default sets follows directly from the 
Propositions El IH and [HI The two series of default sets of the algorithm both converge in finitely 
many iteration steps to D* if there is no firm in the financial system that is borderline. Lemma 
IA5I in the Appendix shows that the probability for borderline firms in R* is zero from which 
almost sure convergence follows. □ 


By its nature, the Sandwich Algorithm converges to D* from both directions which doubles the 
computation and makes the algorithm somewhat inefficient from a computational point of view. 
On the other hand, the algorithm computes an exact solution in finitely many iteration steps 
without wasting time on “Trial-and-Error”. In contrast to the Trial-and-Error Algorithms, the 
drawback of the Sandwich Algorithm is that the convergence of the procedure cannot be ensured 
when borderline firms are present in the system. To overcome this problem, we recommend for 
practical purposes to apply the idea of a lag value in the Sandwich Algorithm as well. 

ALGORITHM 10 (Modihed Sandwich Algorithm). Set I > 2. 

1. Determine ^ and R*^ as well as their corresponding default sets L)(^) and DfR^). 
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_ ^ , 

2. For k>l, calculate the iterates R and R* using one of the Algorithms{^ 0 or0 and the 

_ ^ , 

corresponding default sets DfR ) and DfRf). 

3. If DfR^) = DfR!^), stop the algorithm, set D* = D(R^) and calculate the pseudo solution 
that belongs to D* following Definition^ Else, if k > I and 

|L»(R^)| - \D(r'")\ = ... = |Z?(R^-'+^)| - |L>(R*'"'+^)|, (99) 

_^ 

calculate the pseudo solution belonging to DfR ) and stop the algorithm if it solves the 
Equations ([5]) and ([6]) . Else, set k = k + 1 and go back to step [H 

_ ^ , 

The modification consists of interrupting the algorithm if the default sets D(R ) and DfRf) 
for both iteration directions are not identical but stay constant for I consecutive times. If I is 
chosen large enough (e.g. I > 5) and (IM|) holds, this is a strong indication that at least one 
firm in the system is borderline and that the convergence of both series is not given. In this 
situation, a check whether the default set has already been reached is suitable. 


5. Simulation Study 

Aim of this section is to confirm the theoretical findings in the former sections by simulation. 

In particular, we focus our considerations on the following subjects: 

(i) Investigate the trade-off when choosing a lag value I in the Trial-and-Error Algorithms of 
Section KT[ The result of this part contains ‘optimal’ lag values for each algorithm and 
will be used in the next part. 

(ii) Investigate the algorithm efficiency of all the presented algorithms in the Sections [3] and [H 
For every different technique to calculate the next iterate (Picard, Elsinger, Hybrid), we 
compare the three types of algorithms (non-finite, Trial-and-Error, Sandwich) with each 
other. 

Before presenting the results of our study, the used financial systems are further specihed. 


5.1. General Structure of the Financial Systems 

For the system size n we chose six different values, viz. n G {5,10,25,50,100,200}. A system 
with only 5 or 10 firms can be considered as relatively small whereas networks with n = 25 or n = 


50 are r egarded as medium-s i zed. S mall systems are investigate d for example in iGourieroux et al 


( 2012l i. Rogers and Veraart ( 2013l i and Elsinger et al. ( 2006a) where the size w as 5 , 6 and 10 


firms respectively. Examples of medium-sized systems are lAcemoglu et al.l (l2013l i aud iNier et aJ 
( 2007l i where the size was 20 and 25 firms respectively. Further, we added networks with 100 


and 200 firm s into o ur study to also in clud e larger systems. Existing studies for such sizes are 
Elliott et al. ( 2013l h Cont et al. ( 2O10l ) and Gai et al. (2011) that entailed networks with 100, 
125 and 250 firms, respectively. 

There are some empirical studies dElsinger et al.l (l2006bl h iGai and KaDli^ (l2O10l H that in¬ 
vestigated system sizes of about n = 1000. We believe that for practical purposes such large 
systems are not of interest which is why we did not take values n > 200 into account. However, 
it is expectable that the results obtained for our system sizes also hold for networks with more 
than 200 firms. 
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The next input parameters to define are the asset and debt values. To keep it simple, we 
assumed in every simulation scenario of this study for the exogenous assets a value of 1 for each 
firm, i.e. a = (1,... , 1)* € M"'. For firm z’s nominal debt value, we set a fixed value di = d for 
all i £ Af and added a random variation to each debt value in order to get differing setups which 
leads to 

d = (di,..., dnf + (ei,..., SnY £ K”, (100) 

where the e* are independently normally distributed with mean value 0 and standard deviation 
0.5, i.e. Ei ~ A^(0, 0.25). Note that in case of shocks with e* < di we set d* = 0 to avoid negative 
liabilities. 

When constructing an ownership matrix, the degree of ownership that is operationalized by 
the expression integration can provide some crucial information. 

DEFINITION 5. Consider a financial System T = (a, d, M®). The debt integration 
level is defined as the maximum eolumn sum ofM^, i.e. 

n 

= max V mY = IIM'^II. (101) 

i=i 

Analogously, 

n 

i/® = max Vm®. = ||M®|| (102) 

ieAf ■' 

is ealled the equity integration level. 


The integration level is hence a measure of the extent of cross-ownership in eithe r the debt or 


the eq uity component and its definition is based on the one given in the work of lElliott et al 
( 2013l l. Because of Assumption [2l it follows directly that £ [0,1). 

The integration levels and obviously do not specify the single entries of the ownership 
matrices. For this purpose we will limit our consideration in the following to somewhat regular 
structures of the matrices. 


DEFINITION 6. An ownership matrix M is called 

• a ring ownership matrix if in every eolumn only one entry is larger than 0 and 

• a complete ownership matrix if every entry, except for the diagonal entry, is larger than 
0 and of the same size. 

Further, let M 6e a ring ownership matrix and M 6e a complete ownership matrix. A A-convex 
combination of M and M is defined as the matrix M with entries 

Mij = XMij -|- (1 — A £ [0,1]. (103) 


The concepts of ring and complete matrices and of convex combinations are originally used in 


^ ' 

Acemoglu et al.l ( 20131 1. If is a ring matrix this means that every firm has only one creditor 


within the system, and only one shareholder if M® is a ring matrix. Without loss of generality 
we assume that firm i -|- 1 is the creditor (shareholder) of firm i for i = 1,..., n — 1 and that 
firm n is the creditor (shareholder) of firm 1. When (M®) is a complete ownership matrix, 
debt (share) proportions are equally distributed between the re — 1 firms. The lower A is chosen, 
the more equal are the entries of the corresponding convex combination. 
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EXAMPLE 1. For a system of size n = 4:jwe assume a debt integration level oj^u^ = 0.9 
and set A = 0.5. The ring ownership matrix the complete ownership matrix and the 
X-convex combination matrix then are 




/O 

0 

0 

.9\ 


/o 

.3 

.3 

.3\ 


/o 

.15 

.15 

.6\ 

.9 

0 

0 

0 

iv/fd_ 

.3 

0 

.3 

.3 

nnr! A/fd _ 

.6 

0 

.15 

.15 

0 

.9 

0 

0 

, iVl — 

.3 

.3 

0 

.3 

CLTld IVl — 

.15 

.6 

0 

.15 

VO 

0 

.9 

0/ 


V3 

.3 

.3 



V.15 

.15 

.6 

0/ 


(104) 


Consider two financial systems T = (a, d, M® ) and iF = (a, d, M®) with correspond¬ 
ing integration levels and Due to the regular structure of the ownership matrices 

we can say that a system T is more debt-integrated than iF, if and only if In the same 

way, we define that a system is more equity-integrated. 

With the former definitions, the following parameters are needed for the simulation of a 
financial system: n, d, and A, where we will use the same A to define the debt and 

the equity ownership matrix according to ()103l) . A simulated system is the financial system 
J- = (a, d, M'^, M®), where the parameters n, d, and A are used to define a, and M® 

and where the liabilities d are a realization of the random variable in (jlGOp . 


5.2. Effect of the Lag Value 

As mentioned in Section HU the smaller the lag value I is chosen in the Trial-and-Error Algo¬ 
rithms, the higher is the chance that the first possible default set is not the actual D*. This 
results in unnecessary computation steps to reach the real default set. On the other hand, if I 
is taken as very high, say / = 5 or higher, there are many iteration steps in the algorithm that 
are possibly not needed. For this reason, we wanted to investigate this trade-off situation by 
determining the error rate for the first potential default set. 

Assume that for a given parameters n, d, i^® and A, we have generated N simulated 
systems. For every system we determine for the Trial-and-Error Picard (TP), the Trial-and- 
Error Elsinger (TE) and the Trial-and-Error Hybrid Algorithm (TH) for a lag value I > 2 
the hrst potential default set D^p{l), Z):^g(/) and D^.^{1) where j = 1,N. In case of the 
Trial-and-Error Picard Algorithm we define 



1 , 

0, else. 


(105) 


and analogously £^rj.^{l) and £^^^(0 for the TE and TH Algorithm, respectively. The error rate 
for the TP Algorithm for the lag value I is then given by 

1 ^ . 

eTp(0 = ^ X]^Tp(0 € [0,1]. (106) 

i=i 

In the same way the error rates ste{1) and eTH(0 are dehned. 

For the investigation of the error rate we chose di = d = 1.5 in (llOOl) as the debt value. The 
debt integration values where € {0.9,0.5,0.1}, which we considered as systems with high, 
moderate and low debt cross-ownership. Similarly, we took p® € {0.45,0.25,0.05} for equity 
integration, where each value is half the associated debt integration. The justification for this 
approach is that equity cross-ownership is probably commonly less pronounced than debt cross¬ 
ownership. Further, we wanted to avoid possible cross-ownership entries lager that 0.5 since this 
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would mean that a firm is owned by majority by another firm in the system. Each equity and 
debt integration value was combined with each other which results in 9 possible system settings. 
Beyond that, all 9 settings were investigated three times, where the structure parameter A took 
the three possible values A G {0,0.5,1}, i.e. systems with only ring ownership matrices, systems 
with only complete matrices and systems with a 0.5-convex combination were considered. In 
total, the combination of the parameters and A leads to 27 different settings and for 

every setting N = 1000 simulated systems were generated. The error rates were calculated for 
the three algorithms for the lag values I = 2,... ,7. Repetitive simulations with N = 1000 
showed that the error rates are fairly stable for different simulation runs which is why we viewed 
the number of 1000 repetitions as reliable. We used the Decreasing Trial-and-Error Algorithm 
defined in Algorithm [3 but simulations with the Increasing Trial-and-Error Algorithm showed 
very similar results. 

Table 1: Error rates ^Tp(O) ^te( 0 in percentage points for the Decreasing TP, TE 

and TH Algorithm for I = 2,... ,7. Mean values over all three values of A and all values 
of and p® are shown for each combination. The last three rows of the table show 


overall mean 

error 

rates over 

all considered system 

sizes. 




1 = 2 

1 = 3 

1 = 4 

1 = 5 

1 = 6 

1 = 7 


TP 

9.652 

3.111 

1.170 

0.500 

0.230 

0.104 

n = 5 

TE 

3.452 

0.718 

0.159 

0.037 

0.022 

0.007 


TH 

0.156 

0.015 

0.004 

0 

0 

0 


TP 

4.574 

1.452 

0.529 

0.196 

0.089 

0.048 

n = 10 

TE 

1.256 

0.111 

0.015 

0.004 

0 

0 


TH 

0.233 

0.004 

0 

0 

0 

0 


TP 

8.608 

2.730 

1.063 

0.455 

0.207 

0.067 

n = 25 

TE 

2.819 

0.348 

0.041 

0.004 

0 

0 


TH 

0.533 

0 

0 

0 

0 

0 


TP 

10.596 

3.122 

1.111 

0.482 

0.189 

0.078 

n = 50 

TE 

3.300 

0.389 

0.041 

0.007 

0 

0 


TH 

0.537 

0.011 

0 

0 

0 

0 


TP 

11.359 

3.385 

1.211 

0.478 

0.2 

0.093 

n = 100 

TE 

3.233 

0.411 

0.044 

0.007 

0 

0 


TH 

0.404 

0.007 

0 

0 

0 

0 


TP 

11.485 

3.230 

1.111 

0.419 

0.167 

0.063 

n = 200 

TE 

2.870 

0.322 

0.022 

0.007 

0 

0 


TH 

0.267 

0.007 

0 

0 

0 

0 


TP 

9.379 

2.838 

1.033 

0.422 

0.180 

0.075 

all n 

TE 

2.822 

0.383 

0.054 

0.011 

0.004 

0.001 


TH 

0.355 

0.007 

0.001 

0 

0 

0 


The hrst observation was that the structure of the ownership matrix, that is the choice of A, 
had no severe effect on the error rates. This results in comparing the error rates for systems 
with given n, and for the three choices of A. In case of n = 100 where moderate debt and 
low equity cross-ownership was present = 0.5, = 0.05) and a lag value I = 1, the largest 
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absolute difference is documented for the Trial-and-Error Picard Algorithm between complete 
and ring ownership matrices (9.6% to 3.2%). In the large majority of possible combinations, the 
difference was much smaller which is why we concluded that the ownership structure itself does 
not affect the error rate in an essential way. For this reason, we summarized the three values of 
A and calculated the mean error rates for every combination of n, and over all A for the 
further results. 

In Table [U the results of the simulation to investigate the effect of the system size are snm- 
marized. The error rates are calculated as the mean over all possible combinations for debt and 
equity integration for every value of n. We observe that the network size only slightly affects 
the error rates, since for the same lag value they are relatively close for all system sizes. The 
only exception are systems with n = 10 firms, where the error rate is smaller compared to the 
other systems. In the last three rows of the table, an overall impression of the error rates shows 
that even for the TP Algorithm and a lag value of / = 2, the error rate is not higher than 10% 
in total. We also detected that the error rates for the TE Algorithm are much smaller and even 
more so for the TH Algorithm. For increasing lag values, the error rates quickly diminish in size 
for all considered methods. 


Table 2: Error rates e^p(Z), e^p(Z) and in percentage points for the TP, TE and TH 

Algorithm for I = 2,3. For the debt integration level low, moderate and high 
integration is defined as 0.1, 0.5 and 0.9, respectively. For the equity integration level 
zz® the corresponding levels are defined as 0.05, 0.25 and 0.45. Mean values over all three 
values of A and all values of n are given for each combination. Further, overall mean 
error rates for within each debt and equity integration level are shown in an additional 
column and row. 





1 = 2 





1 = 3 






z^® 





z^® 





low 

mod. 

high 



low 

mod. 

high 


TP 

low 

2.60 

4.65 

7.64 

4.96 

low 

0.22 

0.54 

1.41 

0.72 


mod. 

8.53 

19.64 

35.83 

21.33 

mod. 

1.72 

5.58 

14.62 

7.31 


high 

1.31 

2.47 

1.75 

1.84 

high 

0.14 

0.73 

0.57 

0.48 



4.11 

8.73 

15.14 

9.38 


0.69 

2.29 

5.53 

2.84 









z^® 





low 

mod. 

high 



low 

mod. 

high 


TE 

low 

2.16 

2.34 

2.51 

2.34 

low 

0.19 

0.19 

0.18 

0.19 


mod. 

6.24 

6.07 

5.60 

5.97 

mod. 

1.07 

0.91 

0.86 

0.95 


high 

0.29 

0.13 

0.06 

0.16 

high 

0.03 

0.02 

0 

0.02 



2.90 

2.84 

2.72 

2.82 


0.43 

0.37 

0.35 

0.38 




z/® 



ly^ 


z/® 





low 

mod. 

high 



low 

mod. 

high 


TH 

low 

0.02 

0.13 

0.24 

0.13 

low 

0 

0 

0 

0 


mod. 

0.22 

0.96 

1.59 

0.92 

mod. 

0 

0.01 

0.04 

0.02 


high 

0.04 

0 

0 

0.01 

high 

0 

0 

0 

0 



0.09 

0.36 

0.61 

0.35 


0 

0.01 

0.01 

0.01 
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To assess the influence of the debt and the equity integration level, the mean error rate over 
all considered system sizes are taken and listed for each possible combination of the integration 
level, as shown in Table [5J If the debt integration level increases from low to moderate, the error 
rates increase as well. This can, for example, be seen when comparing the mean error rates for 
the debt integration levels 0.1 and 0.5 over all equity integration levels in the last column in 
Table [2] for every lag value. For the Trial-and-Error Picard Algorithm, the error rate increases 
from 4.96% to 21.33%; for the other algorithms we observe similar results. A further increase of 
the debt integration from 0.5 to 0.9, however, has the reverse effect since the error rates decrease 
in this case. Again we take the TP Algorithm as an example where the error rates shrinks from 
21.33% to 1.84%. A possible explanation for this behavior could be that for the TP Algorithm 
the number of needed iteration steps to converge to the solution was always highest for the 
combination = 0.5 and = 0.45 (data not shown). Hence, the convergence speed in these 
situations is very slow which explains why the first potential default set is often not the actual 
default set. When the mean error rates for every equity integration level averaged over all values 
of is examined, we observed that, except for the TE Algorithm, the error rates increase for 
increasing integration levels. In case of the TP Algorithm, we have an increase from 4.11% to 
8.73% to 15.14% for = 0.05, = 0.25 and = 0.45, respectively. The error rates for the 

TE Algorithm stay approximately constant (2.90%, 2.84%, 2.72%). 

Our overall conclusion of this part of the simulation study is that the choice of the algorithm 
and the lag value I has the strongest effect on the error rate, i.e the error rates quickly decrease 
for greater lag values. The TP Algorithm has, as expected, the highest overall error rates, much 
higher than the TE and the TH Algorithm. What also affects the error rate is the integration 
level of the ownership matrices. Eor our simulation setting, it was the combination of moderate 
debt and high equity integration that yielded in the highest rates. The structure of the ownership 
matrices, on the other hand, had no influence on the error rate. Taking the overall mean error 
rates as the main reference, we can state that for the TE and the TH Algorithm a lag value of 
2 is appropriate, since the corresponding error rates are with 2.82% and 0.36% very small. For 
the TP Algorithm, however, for I = 2 we get an overall error rate of 9.38% which is why a lag 
value of 3 with an overall error rate of 2.84% seems more convenient for this procedure. 


5.3. Comparison of Algorithm Efficiency 


Searching for the most efficient algorithm, the main issue is to minimize the calculation effort to 
find the solution R*. In every iteration step of the algorithm, different kind of calculations are 
carried out for the different algorithms. We quantify the calculation costs with the Landau sym¬ 
bol (Big O notation), where for example 0(n) means that the time T{n) to compute a problem 
of size n grows at the rate n. We distinguish between two different types of calculations in our 
considerations. For the first type, a mapping is applied to a given vector. This mapping can ei¬ 
ther be the mapping 4> in ([71) or the mappings and <I>® in (12611 and (14011 . respectively, whereas 
the second type, matrix multiplications are done. In all cases, the most expensive calculations 
are matrix multiplications, whereas the type embodies the solution of a linear equation system 
such as the ones defined in the Algorithms I2AI andl4Al The computa t ional costs of both types 
are between the range of 0{'n?) and 0(n^) (cf. Dahlauist and Biorck ( 20081 ) L 

Keeping the functioning of the algorithms in mind, the Elsinger and the Hybrid Algorithm 
seem to be less efficient than the Picard Algorithm, since in the latter one, no linear equation 
system has to be solved which results in smaller computational costs. However, as we have 
seen in the Propositions [5] and O in terms of iteration steps, the Hybrid Algorithm converges 
faster to the solution compared to the Elsinger Algorithm, which in turn converges faster to 
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the solution than the Picard Iteration. Therefore, a typical trade-off-situation is given between 
computational costs and convergence speed of an algorithm. Note that for a sequence that 
converges to a fixed point R* the convergence rate is called linear, if there exists a c € (0,1) 
such that 


||R* _Rfc-rl|| < ^||R* _Rfc|| 

for all k>0. Since ll^fR*) - ^>(R^)|| < /”’“||R* - R*^|| with = max{||M‘^||, ||M®||} (see 
Lemma 4.1 in Fischer ( 2014l il. linear convergence holds for the Picard Algorithm if instead of 
Assumption [1] the stronger assumption of matrix norms being strictly smaller than 1 is made 
(which means that of all debt and equity a non-zero share is held by a system outsider). The 
properties of the Elsinger and the Hybrid Algorithms, however, made it impossible to prove 
linear convergence (or an even higher convergence rate). The next problem is that the total 
computational cost for one the algorithms is not determinable in general. For these reasons, the 
comparison of the different calculation techniques (Picard, Elsinger, Hybrid) on an analytical 
base seems impossible. 

This is why we measured the time that was needed to execute an algorithm and considered 
this value as the primary outcome of our simulation. Though this measure strongly depends on 
the processor speed and memory capacity of the computer, it allows an objective comparison of 
the different algorithms. The simul ations were conducte d on a computer with 3.2 GHz and 4 
GB RAM, the software used was R ( R. Gore Team ( 20141 11. 

The parameters dehning the financial systems are given in the following. Unlike to the simula¬ 
tion in Section [521 where a fixed debt value d was used, we varied between hve possible the debt 
values and chose d G {1,1.5, 2,2.5, 3}. The set of equity and debt integration levels was extended 
to G {0.025,0.1,0.175,0.25,0.325,0.4,0.475} and G {0.05,0.2,0.35,0.5,0.65,0.8,0.95}, 
hence seven possible integration level respectively. A result in Section 15. 21 was that the structure 
of the ownership matrix does not influence the error rates, which is why we only took complete 
ownership matrices into account for this simulation, i.e. A = 0 in (11031) . Together with the six 
considered systems sizes (n G {5,10,25,50,100,200}), this new setting results in 6-h-7-7 = 1470 
different settings. Again, for each setting, N = 1000 simulated systems were generated for every 
parameter combination. Eor every simulated systems we applied all 15 algorithms presented in 
the former section and documented the runtime for every procedure to find the solution R*. 
We used the Picard, the Elsinger and the Hybrid Algorithm (Algorithms [H [3| and [6]) with both 
versions, i.e. the decreasing and the increasing version. Eurther, the Trial-and-Error versions 
were considered, again both the decreasing (Algorithm [7]) and the increasing version (Algorithm 
[8]) with lag values of / = 3 for the Picard versions and I = 2 for the Elsinger and the Hybrid 
versions. The choice of the lag value is a result of the simulations in Section 15.21 Note that 
minimizing the error rate to an appropriate value is not necessarily equivalent to minimizing 
the runtime of the algorithms. For these reasons, we compared for every considered scenario 
the runtime of the Trial-and-Error Algorithms using lag values from I = 2 to I = 5. The results 
(not shown here) are that for the Elsinger and the Hybrid versions of the algorithms, the choice 
of f = 2 does not only keep the error rate on a very low level, but also minimizes the runtime. 
Eor the Trial-and-Error Picard Algorithm, the simulation showed that the runtime is almost 
identical for I = 2 and I = 3. However, there is no clear tendency between those choices of 1: 
for some parameter combinations I = 2 lead to smaller runtimes and for some situations this 
was the case for I = 3. Due to this indifference for the Trial-and-Error Picard Algorithm, we 
set I = 3 for these procedures in accordance with the findings of Section [5.21 At last, the three 
versions of the Sandwich Algorithm (Algorithmic]) were also taken into account. The tolerance 
level in all algorithms was set to e = 10“^. 
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Table 3: Mean runtime in seconds for every algorithm over all debt and equity integration values 
and ly'^) and all debt values (d) grouped by system size, algorithm and iteration 
type. For each of the three iteration types, the average runtime of the corresponding 
increasing and decreasing version was calculated, except for the Sandwich Algorithms. 


Algorithm Type 

Iteration Type 



system size n 


5 

10 

25 

50 

100 

200 


Picard 

1.81 

1.90 

2.27 

3.01 

6.19 

21.16 

Non-finite 

Elsinger 

1.86 

2.08 

2.90 

5.71 

24.38 

175.40 


Hybrid 

1.65 

1.84 

2.58 

5.24 

23.06 

164.03 



1.78 

1.94 

2.58 

4.65 

17.87 

120.20 


Picard 

1.29 

1.30 

1.64 

2.81 

8.77 

45.21 

Trial-and-Error 

Elsinger 

1.32 

1.39 

1.91 

3.95 

17.14 

119.38 


Hybrid 

1.57 

1.67 

2.24 

4.49 

19.04 

130.88 



1.39 

1.45 

1.93 

3.75 

14.98 

98.49 


Picard 

1.16 

1.28 

1.79 

3.31 

10.73 

55.08 

Sandwich 

Elsinger 

1.33 

1.54 

2.41 

5.53 

26.18 

192.65 


Hybrid 

1.49 

1.70 

2.55 

5.57 

25.44 

182.91 



1.32 

1.51 

2.25 

4.80 

20.78 

143.55 


Picard 

1.47 

1.54 

1.92 

2.99 

8.13 

37.56 

Overall 

Elsinger 

1.54 

1.69 

2.41 

4.97 

21.84 

156.44 


Hybrid 

1.59 

1.74 

2.44 

5.01 

21.93 

154.55 



1.53 

1.66 

2.26 

4.32 

17.30 

116.19 


An important topic is of course the comparison of new developed techniques (Trial-and-Error, 
Sandwich) with the existing procedure. In Tabled the mean runtimes are listed grouped by the 
size of the financial system as well as the algorithm and the iteration type. We summarized the 
decreasing and the increasing version of every algorithm respectively by calculating the mean 
runtime of both procedures. The runtimes of the decreasing versions were in most situations 
smaller than their counterparts. Ignoring the random structure of d for an instance and calculat¬ 
ing the fixed point, it could be seen that in about 60 % of all considered scenarios, no firm was in 
default which explains the slight ‘overperformance’ of the decreasing algorithms. If we compare 
the mean runtimes over all iteration types, we find that for n = 5 the Sandwich Algorithms 
have the best performance (1.32 s) compared to the Trial-and-Error and the non-finite methods 
(1.39 s and 1.78 s, respectively). For n > 5, the fastest runtimes averaged over the iteration 
types are achieved for the Trial-and-Error procedures. Comparing the different iteration tech¬ 
niques with each other, we find that using the Picard Iteration technique results in a minimal 
computational effort for all considered system sizes. To be more specific, for small financial 
systems (n = 5,10), the Sandwich Picard Algorithm shows the smallest runtime (1.16 s and 
1.28 s, respectively), whereas for medium-sized systems (n = 25,50), the Trial-and-Error Picard 
Algorithm performs best compared to the other algorithm types (1.64 s and 2.81 s, respectively). 
In case of large financial systems, i.e. n = 100,200, the Picard Iteration in its non-finite form 
yields to lowest runtimes (6.19 s and 21.16 s, respectively). In general, it is clearly visible that 
Picard-typed algorithms have the best performance within every Algorithm type. The only 
exception of this trend can be found in the class of non-finite algorithm, where for n = 5 and 
n = 10 the Hybrid Algorithm showed a slightly lower runtime than the Picard Algorithm. In 
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all other situations, however, the Picard Algorithm is superior to the other algorithms. 

Beside the size n, we also investigated the influence of the other parameters that define 
the form of the financial system on the runtime. We observed that increasing debt values d 
result in an increasing calculation effort, see Table 0] in the Appendix for a detailed overview. 
An exception of this tendency represents the Hybrid Algorithm, where the runtime for large 
d begins to decrease again, no matter which algorithm type is considered. The reason for this 
behavior is that the runtime for the increasing versions of the Hybrid Algorithm becomes smaller 
for large n and so does the average of the increasing and the decreasing version of the algorithm, 
that is shown in Table 01 A possible explanation is that the Increasing Hybrid Algorithm uses a 
Picard-type technique to determine the next debt iterate (cf. Algorithm 15All . As shown above, 
the Picard iteration results, in particular for large n, in much better runtime performances. 

If the debt integration level increases, we first observe a similar effect on the runtime as for 
the debt values, i.e. the higher the integration level, the higher the runtime. However, this 
monotonicity holds only up to = 0.5 or = 0.65 in most cases. For larger debt integration 
levels, the computational effort decreases again. The reason for this behavior could be that for 
small values of it is very likely that many firms in the financial system are in default. In such 
situations, we observe that only few iteration steps are needed until R* is reached. If the debt 
integration level is very high, the same effect establishes with the difference that many firms 
in the system are solvent. For medium debt integration levels this clear distinction for a firm 
between solvent and defanlt disappears. A consequence is a higher number of needed iteration 
steps which also influences the runtime. Moreover, this interpretation is underlined by the fact 
that for small (firms more likely in default), the increasing versions of the algorithms have a 
better performances, whereas this relationship inverts for large integration levels where the firms 
are more likely to be solvent. For increasing equity integration levels this effect is not visible, 
the runtime increases if increases (results not shown). Since < 0.45 the equity integration 
seems to have a less strong effect on the status of the firms in the system and therefore, the 
effect seen in the debt integration levels does probably not appear. 


6. Summary 


In this article, we gave a snrvey of the existing algorithms (“Picard” and “Elsinger”) for the 
computation of equilibrium prices in a financial syste m in which cros s-hol dings of equity and 
debt a re present. Moreover, we showed how the ideas of Elsinger ( 20091 1 and Eisenberg and Noe 
(1200 il l can be combined to get an iteration procedure (“Hybrid Algorithm”) that is in every 
iteration step closer to the solution R* than the “Picard” and “Elsinger” algorithms. We 
developed new iteration methods based on the information of defaulting and solvent firms under 
a current payment vector. A consequence of these default set-based methods is that the exact 
solution of the system is reached in a finite number of steps, which could not be ensured for 
the existing iteration procedures. Using this new approach yields to two different concepts, that 
we called “Sandwich” algorithms and “Trial-and-Error” algorithms. While for the former type, 
a clear stopping criteria can be defined (at least almost surely), the latter algorithms have the 
drawback that every potential solntion has to be checked for validity. In a simulation study, 
we showed that choosing an appropriate lag value I, the computational effort can be kept to a 
minimum. 

Another simulation showed that essentially less iteration steps have to be performed when 
using the new default-set based techniques. However, the faster convergence concerning the 
number of iterations has its price: In the new algorithms other than the Picard type, potentially 
several linear equation systems have to be solved in every iteration step. This leads to a higher 
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calculation effort for those methods and a result of the empirical investigation of the runtime 
for all algorithms was that in particular for large financial systems the computational costs then 
become higher than for algorithms where no linear equation systems have to be solved. Another 
result of the runtime analysis is that the most efficient iteration technique is of Picard type. 
In the majority of the considered settings those iteration techniques performed best, no matter 
which algorithm type (Non-hnite, Trial-and-Error, Sandwich) was used. One of the main results 
is that the choice of the most efficient algorithm strongly depends on the size n of the hnancial 
system. We observed that for small systems (n = 5,10) the Sandwich Picard technique, for 
medium-sized systems (n = 25, 50) the Trial-and-Error Picard technique and for large systems 
(n = 100, 200) the simple non-hnite Picard technique achieves the best results with regard to the 
minimization of the runtime. Regarding the choice of the tolerance level e, smaller values than 
the used one of e = 10“^ will strongly affect the results of the non-hnite iteration techniques, 
since additional simulations (results not listed here) revealed that an increase of e will lead to a 
disproportionally strong increase of the needed iteration steps and therefore the runtime. One 
consequence could be that for large systems the non-hnite Picard iteration would not be optimal 
anymore, since the hnite algorithm techniques do not depend on e. We are aware of this effect, 
but we think that, for practical purposes, a tolerance level of e = 10“^ is sufficiently small. 

The simulation in Section [5.31 contained only complete ownership matrices (A = 0 in (I103p i. 
It is of potential interest, whether in case of ring ownership matrices or A-convex combinations 
the results lead to the same conclusions. No matter which value of A is chosen, the entries of 
and M® still are uniquely determinable. A potential extension of this assumption would be 
to allow random own ership matrices based on a random network matrix as used for example 
( 2ni3l i. Besides these questions, the main focus for further research should be 


m 
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on generali z ing th e algorithms for systems with more than one seniority level for the liabilities . 
In Fischer (2014), this was done for the non-hnite Picard Algorithm and in Elsingei ( 2009l b 
an extension of the non-hnite Elsinger Algorithm is discussed as well. It would be of interest 
whether and - if yes - how the Hybrid Algorithm can be generalized for a model that allows for 
a seniority structure of debt. 


A. Appendix 


A.l. Proofs and Auxiliary Results 


LEMMA Al. Let || • || be a not necessarily strictly convex norm on M”', and let ^ be a map 
on a nonempty convex and compact set C C M" which is non-expansive with respect to the 
norm-induced metric. The set of fixed points of ^ in C is then nonempty, closed, and either a 
singleton, or uncountable. 


Proof. Much-rehned versions of this result are known (e.g. Bruck ( 1973I H. For convenience, a 
short proof is given. Non-expansiveness implies that is (1-Lipschitz) continuous. The set 
of hxe d points is hence closed, and the Brouwer-Schauder Fixed Point Theorem (e.g. iRudin 
(|l99ll H provides the existence of at least one hxed point. Assume now that x, y G C are two 
distinct hxed points of 4>. For v G C and e > 0, Be(v) = {w G C* : ||w — v|| < e} is a non-empty, 
convex and compact subset of C. For A G (0,1), the intersection 


= -SA||y-x||(x) n R(i_A)||y_x||(y) (108) 

is non-empty (as it contains (1 — A)x -|- Ay)), convex and compact, and it contains neither x, 
nor y. By the triangle inequality, fl Cx^ = 0 for Ai A 2 . Non-expansiveness implies that 
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‘^(Ca) C Cx- By Brouwer-Schauder, there exists a fixed point of in Cx- Hence there exist 
uncountably many fixed points of ‘h in C. □ 


LEMMA A2. Let M G be an ownership matrix that has the Elsinger Property. Then 

/ 9 (M) < 1, where 

/ 9 (M) = max{|Aj| : A* eigenvalue o/M} (109) 


is the spectral radius o/M. 

Proof. A well known result (cf. R.udinI ( 199 il l 1 is that / 9 (M) < ||M|| < 1. In case of ||M|| < 1 
there is nothing to show, so we assume that ||M|| = 1 which is no contradiction to the Elsinger 
Property of M. We will show the claim by contradiction. To this end, assume that /j(M) = 1. 
For the corresponding eigenvalue v is must hold that v 7 ^ 0 and Mv = p(M)v = v. We can 
formulate this equation alternatively as 


(I„ - M)v = On. (110) 

Since M has the Elsinger Property, it follows by Elsinger ( 2009l l. Lemma 1, that (1^ — M) is 
invertible. But that means that there exists no vector v 7 ^ 0 such that (| 110|1 is true. Hence, 
V = 0 which is a contradiction and from which follows that p(M) <1. □ 

LEMMA A3. Let M G be an ownership matrix that has the Elsinger property and for 

which /o(M) < 1. Then (I„ — M)“^ exists and can be obtained via the Neumann expansion." 


(I„-M)-i = ^M", 


( 111 ) 


n=0 


where M° = I„. Consequently, the diagonal entries of (I^ — M) ^ are greater than or equal to 
1 and the other entries are all non-negative. 


Proof. See lRudinI (| 199 il l. 


□ 


Proof of THEOREM [p A p roof of Theorem [J i s necessary because related proofs in ISuzuki 
( 2002l b Gourieroux et al. ( 2012 1 a ndlEischer (|2014l l rely on stronger matrix conditions than the 
Elsinger Property, while Elsinger ( 20091 ') considers an equation system which slightly differs from 
(0) and First, note that ([5|) and ([U]) can only have non-negative solutions. This is shown in 
Lemma 3.5 of Fischer ( 2014l l under stricter matrix conditions, but because of Lemma [A3l of this 
paper, it is straightforward to see that the proof works in the same manner under the Elsinger 
Property. The interval [Rsmaih R-great] is convex and compact, and ‘h(R) is continuous in R. By 
Lemma [Hand the Brouwer-Schauder Eixed Point Theorem, it follows that at least one solution 
ex ists. Fu r therm ore, as in Eq. ([7|) is a non-expansive mapping. This follows from Lemma 4.1 
of [Fischer] (j2014l l. where a strict contraction property is shown under stricter matrix conditions, 
but again it is straightforward to see how the corresponding proof implies non-expansiveness 
under the Elsinger Property for all ownership matrices. Since it follows from Proposition [TO] 
that there can be a maximum of 2”" possible solutions of ([5]) and Q, uniqueness follows from 
Lemma [Hand Lemma lATl in the Appendix. □ 


LEMMA A4. Let M G be an ownership matrix as in Lemma 

and the matrix A G be defined as 


(A)., = 


1) if i = j o^nd i G A/q, 
0 , else. 


such that Mq C M, 


( 112 ) 
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Then it holds that 


Proof. Note that 


I„ - AMA)"^A < A(I„ - M)~^A. 


(I„ - A)'= = {In - A) for fc G N 
and that = (I„ — A)® = I„. Using Lemma lA3l we have that 


/ OO 




(In - AMA)-iA = ^(AMA)^ A 

\n=0 / 


= (In + AMA + AMAMA + AMAMAMA + .. .)A 
= A +A(M + MAM + MAMAM + ...)A 

<M2 <M3 


< A + A ( ^ M" ) A 


\n=l 


A ( ^ M” j A 


'^n=0 


= A(In - M)-iA. 


(113) 

(114) 


(115) 


□ 


LEMMA A5. Let the random variable a have a have a density with respect to the Lebesgue 
measure on (M(J")"’. The set of all a for which the system solution contains at least one borderline 
firm, i.e. one i € M such that r* = di and s* = 0, then has measure zero. 

Proof. First, note that it suffices to show the claim for a set A{I) of all a for which r* = di 
and s* = 0 for each z G / C A/", since the number of subsets of {1, is finite and a finite 

union of sets of Lebesgue measure zero has Lebesgue measure zero. We first sh ow that A(/) i s 
a Borel set and hence Lebesgue measurable. For this, note that it is shown in iFisched (j2014l ') 
that the mapping 4* : a i—R.*(a) that maps any price vector of the exogenous assets onto the 
corresponding solution of ([5]) and ([6]) is Borel measurable. Let now H{I) denote the 2(n — |/|)- 
dimensional hyperplane in for which 


H{I) = {(r , sj* G R^” : Vi = d and Sj = 0 for all i G /}. 
Clearly, H{I) is a Borel set. One obtains 

A{I) = ^-\h{i) n (R([)^”), 


(116) 


(117) 


which must be Borel-measurable, too. Observe now that if a 2 ^ ai (a 2 strictly larger than 
ai in all components), then <I>” 2 (R-) > ^ai(^) non-negative R. Hence, by the Picard 

Iter ation, R * (a 2 ) > R*(ai). From Eq. ([5]) and dH), it follows that r -|- s = a -|- M'^r -|- M®s (see 
also Fischer (2014)). Therefore, if a 2 » ai , then r*(a 2 ) -|-s*(a 2 ) > r*(ai) -|-s*(ai) for all i € M, 
which is a contradiction to ai, a 2 G A{I). Thus, since r*(a) -|- s*(a) = di for a G A{I) and i £ I, 
a 2 ^ ai can hold for no pair ai,a 2 G A{I). This means that A{I) bears some resemblance 
to a Pareto set. It follows that the set A{I) intersects any straight line parallel to the vector 
(1,..., 1)* G R” either once, or not at all. As such, and since the Lebesgue measure is rotation 


invariant, the problem reduces now to the one which is shown in Lemma lA61 


□ 
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LEMMA A6. Let Q be a Borel set in R”’ such that < 1 for any ui € R"’ where 
Qu = {x G R : {x,u}^Y G Q}. Then Q has Lebesgue measure zero. 


Proof. Let Xmi'na G N, denote the Lebesgue measure o n R”^. For any Borel set Q, it follows 
from the definition of product measures (e.g. Billingsley ( IQQSli l and from = Ai An_i that 


An(Q) 


Al 1 (w). 


(118) 


Since Xi{Qu)) = 0, the result follows. 


□ 


Table 4; Above: Mean runtime in seconds for every algorithm over all debt and equity integration 
values and and system sizes (n) grouped by the debt values (d). Below: Mean 
runtime in seconds for every algorithm over all equity integration values debt 

values (d) and system sizes (n) grouped by the debt integration values 


Algorithm 



P 

E 

H 

TP 

TE 

TH 

SP 

SE 

SH 

1 

5.43 

31.21 

33.87 

9.51 

22.71 

26.78 

10.94 

30.90 

35.17 

1.5 

5.61 

35.73 

36.59 

9.95 

25.03 

28.61 

11.86 

38.27 

39.72 

d 2 

5.98 

36.46 

34.70 

10.31 

25.00 

27.56 

12.55 

40.92 

39.19 

2.5 

6.46 

36.67 

31.69 

10.53 

24.65 

26.23 

12.97 

41.04 

36.29 

3 

6.81 

36.87 

28.49 

10.56 

23.53 

24.07 

12.80 

40.22 

32.67 

0.05 

4.04 

17.24 

20.85 

8.86 

13.21 

18.79 

9.13 

15.25 

20.49 

0.2 

5.34 

26.56 

27.01 

9.41 

17.04 

23.52 

10.42 

22.60 

28.80 

0.35 

6.38 

36.29 

32.89 

10.14 

22.30 

27.12 

11.95 

32.59 

36.41 

0.5 

7.40 

47.30 

40.33 

11.44 

30.18 

31.53 

13.98 

47.05 

44.93 

0.65 

7.45 

55.41 

46.72 

11.75 

35.70 

34.44 

15.24 

61.52 

51.42 

0.8 

6.21 

39.27 

37.51 

10.23 

28.01 

28.20 

13.31 

50.73 

41.69 

0.95 

5.58 

25.65 

26.16 

9.37 

22.84 

22.95 

11.53 

38.17 

32.53 


A.2. Additional Tables and Simulation Results 

The notation in the tables in this section is as follows. The names of the algorithms in the table 
are composed out of their iteration type (“P” for Picard, “E” for Elsinger and “H” for Hybrid) 
and their direction (“D” for decreasing, “I” for increasing). If the prefix “D” or “I” is omitted, 
the mean value of the corresponding increasing and decreasing version is shown. The additional 
prefix “T” denotes the Trial-and-Error version and “S” denotes the Sandwich version of the 
algorithm. 
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Table 5: Median of the iteration (calculation) steps for every algorithm over all debt and equity 
integration values and v^) and debt values {d) grouped by the system size (n). A 
calculation step is defined as the solution of a linear equation system, which is done for 
example in Algorithm I2AI in every iteration step. An iteration step is defined as the 
step from the fc-th iterate to for A: > 0, no matter which algorithm is used. 





system size n 




5 

10 

25 

50 

100 

200 

DP 

6 

7 

8 

8 

9 

9 

IP 

8 

9 

10 

10 

11 

11 

DE 

4(5) 

5 (6) 

5 (8) 

5 (9) 

5 (10) 

6 (12) 

IE 

4(6) 

5 (7) 

6(9) 

6 (11) 

7(12) 

7(13) 

DH 

2 (5) 

3(7) 

3(8) 

3(9) 

3(11) 

3(11) 

IH 

2 (4) 

3(4) 

3 (4) 

3 (5) 

3(6) 

3(8) 

DTP 

2 (1) 

2 (1) 

2 (1) 

3(1) 

3(1) 

3(1) 

ITP 

2 (1) 

3(1) 

3(1) 

4(1) 

4(1) 

5 (1) 

DTE 

1 (3) 

1 (3) 

1 (5) 

1 (5) 

2 (5) 

2 (7) 

ITE 

1 (3) 

1 (4) 

2 (5) 

2 (7) 

2 (7) 

3(8) 

DTH 

1(4) 

1(4) 

1 (6) 

1 (6) 

2 (7) 

2 (8) 

ITH 

1 (3) 

1(4) 

2(4) 

2 (5) 

2 (6) 

2 (7) 

SP 

1 

1 

2 

2 

3 

3 

SE 

0(4) 

1(4) 

1 (6) 

1 (8) 

2 (8) 

2 (12) 

SH 

0(4) 

1(4) 

1 (5) 

1 (6) 

1 (8) 

1 (7) 
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